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Abstract 

This  report  is  concerned  with  the  theoretical  and  experimental  study 
of  subharmonic  generation  in  acoustic  systems.  The  generalized  formulation 
for  lumped  systems  is  considered  for  the  case  of  three  oscillators  coupled 
through  a  nonlinear  element.  The  resulting  analysis  indicates  that  the  high 
frequency  oscillation  is  unstalle  and  its  energy  can  be  diverted  to  low  fre¬ 
quency  oscillations;  that  is,  subharmonics  are  generated.  Based  on  this 
mathematical  model,  the  conditions  for  sustaining  subharmonics  are  deduced 
in  terms  of  the  system's  parameters.  For  distributed  systems,  the  queai- 
linear  solution  of  the  nonlinear  wave  equation,  obtained  by  a  two-variable 
(two-timing)  perturbation  method,  shows  that  the  nonlinear  property  of 
medium  provides  a  coupling  effect  for  signals  of  different  frequencies. 

In  a  multi-resonant  system,  certain  subharmonic  modes,  determined  by  a 
selection  rule,  can  be  excited  by  an  external  source.  The  present  theory 
predicts  the  threshold  for  subharmonic  generation  as  a  function  of  the 
losses,  nonlinearity,  and  detuning  parameters  of  the  system.  Experimental 
measurements  of  such  phenomena  in  an  interferometer-type  resonator,  operated 
around  1,5  MHZ  in  water,  are  consistent  with  the  analytic  results,  provided 
that  cavitation  in  the  medium  is  carefully  avoided. 
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Subharmor. ic  veneration  in  a  •''r.a^oner.on  which  has  beao 
observed  in  acoustic  systems  as  a  re su 1  t  of  their  response  to 
lr tense  signals.  The  turrose  of  this  thesis  is  to  develop  a 
simnle  m.athemat ical  code]  for  ths  basic  mechanism  of  suhhar- 
rnonic  veneration  and,  undo **  the  conditions  arranged  in  tre 
laboratory,  to  investigate  such  a  phenomenon  experimental ly 
for  creaking  the  adequacy  of  tne  theoretical  analysis, 

We  introduce  th.2  subject  with  a  study  on  the  behavior 
of  lumped  systems  with  rass  and  spring  elements.  This  par¬ 
ticular  problem,  is  formulated  by  considering  three  couni ed 
oscillators  with  a  sircie  nonlinear  elem.ert  whose  energy 
function  is  rrorortior.al  to  the  product  of  tne  displacement 
amplitudes  lr.  these  trree  oscillators.  In  a  conservative 
system,  the  analysis  result?  in  three  constraints  for  Tne 
response  to  the  disturbance.  That  is,  in  addition  to  energy 
conservation,  there  are  relations  for  the  inbalance  in  energy 
exchange  among  tne  osc t1  latn^r ,  and  the  amplitudes  and  phases 
of  the  response.  These  corstraints  car.  ce  used  to  reduce 
such  a  problem  tc  quadrature.  The  phase  diagram  for  de¬ 
scribing  the  system's  resnnnse  indicates  that  the  high  fre¬ 
quency  oscillation  is  unstable  and  its  energy  can  be  easily 
diverting  to  the  low  frequency  oscillation.  This  tendency 
suggests  that,  for  a  noncons^rvativc  system,  subharmor.ic 
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oscillation  can  be  excited  and  sustained  if  the  energy  sup¬ 
plied  can  balance  out  the  energy  losses  in  the  system. 

Based  on  this  model,  derive  the  threshold  for  subharmonic 
oscillation  in  relation  to  the  nonlinear  coupling,  dissi¬ 
pation.  and  detuning  parameter  of  the  system. 

We  then  discuss  subharmonic  generation  in  distributed 
systems  by  examining  the  solutions  of  the  nonlinear  wave 
equation.  We  find  that  the  nonlinear  property  of  the  medium 
car  provide  a  coupling  effect  for  signals  with  distinct  fre¬ 
quencies.  Tn  free  space,  an  intense  acoustic  wave  will  form 
a  shock  wave  and  no  subharmonic  can  be  generated.  But  in  a 
closed  system,  such  as  a  resonator,  there  is  the  possibility 
for  subharmonic  modes  to  be  excited.  The  mathematical  analy*» 
sis  for  mode-coupling  in  a  resonator  shows  some  similarities 
with  the  nonlinear  coupling  of  the  three  oscillators  in  the 
lumped  system.  Since  a  resonator  generally  possesses  many 
modes,  the  analysis  has  been  extended  to  determine  a  selec¬ 
tion  rule  for  which  modes  will  be  excited  at  a  given  condi¬ 
tion. 

Experimental  investigation  of  subbannonic  generation 
is  carried  out  ir.  an  interferometer  type  resonator.  Water 
is  used  as  the  medium  that  supports  the  intense  standing 
wave  fo»*  inducing  the  nonlinear  coupling  effect.  Because 
cavitation  is  carefully  avoided  by  degassing  and  filtering 
the  liquid,  we  are  able  to  obtain  results  that  are  consistent 
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with  our  theoretical  model. 

From  the  results  of  this  research,  we  reach  the  con¬ 
clusion  that  the  subharmon5.c  generation  phenomenon  is  attrib¬ 
uted  to  nonlinearity  in  the  system  and  its  sustenance  also 
depends  on  the  loss  factor  and  detuning  parameter  associated 
with  the  particular  subharmonic  modes  to  be  excited.  In  a 
water-f illed  resonator  whose  linear  dimension  is  larger  than 
serval  wavelengths,  v/e  have  found  that,  in  the  absence 
cavitation,  the  generation  of  subharmonics  is  mainly  due  to 
the  nonlinear  property  of  the  medium. 
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Chapter  I 
INTRODUCTION 


The  subharmonic  phenomenon  has  been  studied  in  this 
research,  with  particular  emphasis  on  its  existence  in  acous¬ 
tical  systems.  We  shall  first  discuss  the  nature  of  sub¬ 
harmonics  and  review  some  related  work  reported  in  the  liter¬ 
ature.  We  shall  then  define  the  purpose  of  doing  such  a 
stuay . 

1.  Subharmonics  and  superharmonics 

In  a  nonlinear  system,  a  response  at  other  than  the 
driving  frequency  is  often  observed.  Most  often,  the  fre¬ 
quencies  of  the  response  in  such  a  system  are  related  to  the 
driving  frequency  by  an  interger  multiplier,  and  these  are 
known  as  harmonic  responses.  Eut  under  certain  conditions, 
a  response  with  a  frequency  less  than  the  driving  frequency 
can  also  appear.  To  distinguish  between  these  two  phenomena, 
the  response  with  frequency  less  than  driving  frequency  are 
known  as  subharmonics  and  the  others  as  superharmonics. 

Though  both  subharmonics  and  superharmonics  are  phenom¬ 
ena  caused  by  non! inear it ies  in  a  system,  they  have  a  differ- 
1* 

ent  character  .  Superharraonics  are  a  stable  response  that 


*  Numerical  superscripts  corresoond  to  references  listed  in 
the  Bibliography. 
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will  always  occur  when  there  is  a  driving  source  in  the  non¬ 
linear  system*  Their  frequency  component  is  related  to  the 
order  of  the  nonlinearity  of  the  system.  On  the  other  hand, 
in  order  for  subharmonics  to  be  excited,  a  certain  minimum 
driving  signal  strength  is  required.  The  subharmonic  oscil¬ 
lation  is  due  to  an  unstable  property  of  the  system  and 
their  sustained  excitation  depends  not  only  on  the  intensity 
of  the  driving  source  but  also  on  the  proper  phase  relation 
with  the  driving  source.  Except  for  the  degenerate  case  of 
subharmonics  of  one-half,  subharmonics  ordinarily  exist  in 
pairs  such  tnat  the  sum  of  each  frequency  pair  equals  xhe 
frequency  of  the  driving  signal. 

The  study  of  subharmonics  and  super harmonics  in  non¬ 
linear  mechanical  and  electrical  systems  has  been  carried 

2  13 

out  in  some  detail  by  Kayashi  and  Minorsky  .  Stern's-^ 

book  on  nonlinear  systems  analysis  has  extended  the  non¬ 
linear  phenomena  into  multidimension  space  vectors.  Recent- 

4  *5 

ly  ,  Kronauer  and  Musa  have  done  some  general  analyses  on 
the  synchronization  of  subharmonics  with  the  excitation  for 
weakly  coupled  nonlinear  systems. 

Early  studies  of  sub-  and  superharmonic  phenomena  were 
pursued  for  the  purpose  of  minimizing  undersirable  harmonic 
effects  in  mechanical  systems.  These  analyses  generally  end 
up  predicting  the  conditions  for  which  the  system  will  become 
unstable  and  anharmonic  oscillations  will  begin  to  grow. 
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Later,  owing  to  the  application  of  nonlinear  elements  in 
electrical  circuits,  it  was  found  that  harmonic  generation 
can  be  used  as  a  substitution  for  frequency  multiplication 
when  a  direct  means  for  obtaining  such  a  source  is  not 
available.  Since  the  recent  development  of  high  intensity 
laser  beams,  utilisation  of  harmonic  generation  has  even 
been  extended  to  optics.  The  object  of  subharmonics  studies 
has  also  been  expanded  from  merely  understanding  their  exist¬ 
ence  to  research  on  the  property  of  systems.  From  compari¬ 
sons  with  the  results  of  parametric  amplifier  studies^,  sub¬ 
harmonic  generation  seems  to  have  features  similiar  to  the 

down-conversion  amplifier.  According  to  the  Manley-Rowe 

£ 

relation  ,  high  gains  of  such  conversion  should  be  easy  to 
achieve  through  the  subharmonic  generation  mechanism,  but 
its  unstable  character  and  very  narrow  bandwidth  have  limit¬ 
ed  its  practical  application.  Subharmonic  generation  can 
also  be  employed  in  frequency  dividers  if  the  proper  synchro¬ 
nization  can  be  adequately  maintained.  Hov/ever,  recent 
electronic  digital  techniques  have  provided  a  much  simpler 
arrangement.  This  makes  obsolete  the  idea  of  using  the  sub¬ 
harmonic  generation  mechanism  as  a  alternate  frequency 
soui'ce . 

7 

Tucker  ,  in  his  study  of  nonlinear  properties  in  under¬ 
water  acoustics,  has  discussed  some  possible  applications  of 
subharmonics  in  forming  narrow  low  frequency  acoustic  beams. 


r — . 
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In  the  next  section,  we  shall  review  some  of  the  work  relat¬ 
ed  to  subharmonic  generation  in  the  field  of  acoustics. 

2.  Work  on  Subharmonics  in  Acoustics 

The  earlist  work  related  to  subharmonics  in  acoustics 

O 

is  probably  due  to  Helmholtz  ,  He  described  the  phenomenon 
in  which  one  can  sense  a  fundamental  pure  tone  when  only  two 
of  its  harmonics  are  present.  This  has  been  recognized  later 
as  the  beating  phenomenon  whose  real  cause  is  not  the  inter¬ 
action  of  sound  waves  hr  l  the  r.oni;_ne«:r  response  of 

Q 

the  human  ear.  Lord  Rayleigh'  deduced  solutions  for  a 

3  2 

Duffing  equation  (u  +  ku  +  k*uy  +  n  u  =  0)  and  indicated 
that,  in  the  presence  of  a  driving  force,  a  response  of  less 
than  the  frequency  of  the  driving  source  can  be  excited. 

Work  of  this  kind  has  been  studied  extensively  in  the  field 
of  nonlinear  mechanics  ever  since. 

The  extensive  use  of  sound  for  underwater  communica¬ 
tion,  initiated  in  World  War  II,  has  provided  a  great  deal  of 
research  into  the  generation,  transmission,  and  reception  of 
sound  waves  in  a  liquid.  One  of  the  major  efforts  is  to 
maximize  the  range  of  acoustic  information  transmission. 
However,  there  are  some  limitations  to  such  efforts.  One  of 
them  is  the  cavitation  problem.  When  the  sound  wave  pressure 
reaches  too  great  an  intensity,  foggy  acoustic  streamers 
form  in  front  and  on  the  surface  of  a  transducer  (this  phe- 
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nomenon  is  what  we  now  call  "gaseous"  cavitation).  Such 
streamers  produce  noise  detectable  by  hydrophones.  This 
noise  has  a  line  spectrum  containing  subharmonic s^.  Sub¬ 
harmonics  therefore  become  one  of  the  interesting  subjects 

in  the  study  of  cavitation  phenomena.  Some  investigators^’ 
12 

even  suggest  that  the  subharmonic  can  be  regarded  as  an 
indication  for  the  occurrence  of  cavitation. 

Tucker  has  pointed  out  some  possibi] ities  in  under¬ 
water  acoustics  for  the  utilization  of  the  nonlinear  char¬ 
acter  of  the  med  urn.  One  such  experimental  investigation 
is  the  generation  of  subharmonic s1 ^ .  A  focused  standing 
wave  system  is  employed  during  the  investigation  for  build¬ 
ing  up  a  very  strong  acoustic  field.  When  the  acoustic 
pressure  reaches  a  certain  level,  the  subharmonic  has  been 
observea.  Such  a  threshold  appears  to  be  correlated  with 
the  gas  content  of  the  liquid.  V/hen  boiled  water  is  used, 
those  investigators  reported  that  no  subharmonics  were 
observed. 

14 

Eller  and  Flynn  have  focused  their  study  on  the  non¬ 
linear  problem  of  subharmonics  excited  by  gas  bubble  oscil¬ 
lations  in  the  liquid.  Their  analytical  results  show  that 
a  bubble  can  be  excited  to  its  resonance  when  it  is  subjected 

to  an  acoustic  field  with  twice  the  natural  frequency  of 

1 S 

such  a  bubble,  Neppiras  did  some  experimental  work  along 


this  line  by  injecting  bubbles  of  controlable  size  into  a 
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liquid.  The  results  indicate  that  subharmonics  are  easily 
excited  and  their  intensity  is  strong  when  the  condition  of 
Eller  and  Flynn  described  above  is  established.  They  con¬ 
cluded  that  subharmonics  are  generated  by  bubbles  through 
the  cavitation  process. 

On  a  separate  front,  subharmonics  have  also  been  an 
interesting  subject  under  discussion  in  the  field  of  ultra¬ 
sonics.  While  studying  light  diffraction  patterns  produced 

1.6  17  1 8  IQ 

by  standing  waves  in  a  liquid  *  ‘ ’ “  f  ,  it  is  observed 

that  extra  dots  appear  on  the  regular  diffraction  pattern 

when  the  sound  pressure  reaches  a  certain  level.  Those 

additional  dots  corresoond  to  another  imposed  standing  wave 

whose  wavelength  is  larger  than  that  of  the  original  driv- 

20 

ing  signal.  Cook  has  explained  such  a  phenomenon  by 
assuming  possible  wave  interaction  for  subharmonic  genera¬ 
tion,  but  he  did  not  obtain  a  complete  description  of  the 

relation  between  the  frequency  of  driving  signal  and  the 

21 

subharmonics  observed.  Adler  and  Breazeale  have  inter¬ 
preted  this  phenomenon  differently.  They  concluded  that 
the  vibration  of  the  boundary  of  the  acoustical  interfer¬ 
ometer  parametrically  excites  the  subharmonics  around  the 
order  of  one-half  of  the  original  signal.  However,  such 

a  hypothesis  has  ignored  two  observed  experimental  facts: 

16 

that  subharmonics  of  other  order  have  also  been  observed  ’ 
17  1 fi  IQ 

'*  ’  y  and  that  only  at  a  certain  length  of  the  interfer- 
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16 

ometer  they  can  be  excited  . 

Subharmonic  phenomena  have  been  detected  in  solids 

22  23  24 

when  an  acoustic  wave  is  applied  to  crystals  ’  . 

Luukkala  interprets  these  observations  in  terms  of  a  phonon 
breakdown  hypotnesis.  The  qualitative  phenomenlogical  argu¬ 
ments  he  uses  seem  to  support  the  assumption  that  the  gen¬ 
eral  instability  of  subharmonic  signals  is  caused  by  exces¬ 
sive  energy  going  into  an  additionally  accessible  vibration 
mode.  These  arguments  correctly  predict  the  threshold  for 
one-half  subharmonics  but  are  not  comprehensive  enough  to 

include  three -phonon  interaction, 
oc  26  27 

Dallos  ’  '  reports  observing  subharmonics  in  audi¬ 
tory  systems.  That  such  subharmonics  also  occur  in  pairs 
and  that  sum  of  these  frequency  pairs  equals  the  driving 
frequency  seems  to  indicate  a  similarity  with  other  physi¬ 
cal  systems. 

3.  Scope  of  This  Thesis 

The  purpose  of  this  thesis  is  to  study  the  basic 
principles  underlying  subharmonic  phenomena;  no  attempt 
is  made  to  explain  every  aspect  of  the  existing  observations. 
We  use  a  simplified  mathematical  model  in  hopes  of  obtaining 
some  fundamental  quantitative  relations  among  the  important 
parameters.  In  the  accompanying  experimental  work,  we  not 
only  try  to  confirm  the  theory  but  also  intend  to  clarify 
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some  inconsistencies  among  other  reported  experimental  inves¬ 
tigations.  We  are  concerned  with*  What  is  the  mechanism 
for  subharmonic  generation?  What  is  the  selection  rule  for 
the  appearance  of  a  given  subharmonic  mode?  What  is  a  sys¬ 
tematic  and  practical  way  to  excite  subharmonics? 

In  Chapter  II,  we  study,  in  a  general  sense,  subharmon¬ 
ic  generation  in  lumped  systems.  A  model  of  thi'ee  conserva¬ 
tive  oscillators  coupled  weakly  through  a  nonlinear  element 
is  used  for  the  analysis.  Using  a  first  order  perturbation 
method,  we  have  found  that,  besides  the  energy  conservation 
law,  such  a  system  possesses  two  more  invariants—  enough 
to  reduce  the  problem  to  quadrature.  The  analysis  indicates 
that  the  high  frequency  mode  can,  under  certain  conditions, 
divert  its  energy  to  low  frequency  modes.  By  considering, 
then  the  effect  of  dissipation  in  the  system,  we  have  de¬ 
rived  the  threshold  level  of  excitation  for  sustaining  such 
subharmonic  oscillation.  The  detuning  problem,  that  is  when 
the  frequencies  of  subharmonics  do  not  exactly  match  the 
driving  signal  frequency,  has  also  been  treated. 

We  discuss  subharmonic  generation  in  distributive  sys¬ 
tems  in  Chapter  III.  The  response  of  the  system  has  been 
examined  for  disturbances  with  a  finite  travelline:  soeed. 

A  wave  equation  is  formulated  for  studying  the  effects  of 
all  pertinent  parameters.  We  find  that,  due  to  the  nonlinear 
property  of  the  medium,  there  is  a  coupling  relation  among 
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signals  with  distinct  frequencies,)  In  a  multiresonator,  it 
becomes  mode-coupling  which  has  features  similar  to  the 
coupling  of  oscillators  described  in  Chapter  II,  The  condi¬ 
tion  for  forcing  excitation  of  certain  subharrnonic  modes  is 
then  obtained  analytically  based  on  information  about  the 
boundary  conditions,  detuning  parameter,  and  the  dissipation 
factor  of  the  system. 

The  detailed  experimental  set-up  is  described  in 
Chapter  IV.  The  purpose  of  such  an  investigation  :s  to  seek 
evidence  in  supporting  of  our  hypothesis  about  subharmonic 
generation.  The  liquid  used  in  this  experiment,  water,  has 
been  filtered  and  degassed  in  a  controllable  manner  in  order 
to  get  consistent  results.  Data  are  collected  in  terms  of 
actual  acoustical  quantities  so  that  an  accurate  physical 
interpretation  can  be  made. 

Chanter  V  presents  the  conclusions  of  this  research. 

Y/e  comoare  our  analysis  and  our  experimental  results,  and 
examine  validity  of  our  hypothesis  in  understanding  the 
mechanism  of  subharmonic  generation.  The  direction  for  the 
future  work  and  some  possible  applications  of  the  subharrnonic 
phenomenon  are  suggested. 


{ 
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Chapter  II 

SUBHARMONIC  GENERATION  IN  LUMPED  SYSTEMS 

The  problem  under  consideration  in  this  chapter  is  the 

28 

energy  exchange  among  three  oscillators.  Gilchrist  ,  in 
his  investigation  of  conservative  quasilinear  systems  with 
two  degrees  freedom,  discovered  an  integral  constraint  on 
the  amplitudes  of  the  oscillators.  This  constraint  actually 
is  an  energy  conservation  law  of  first  orde  in  a  pertur¬ 
bation  expansion.  An  additional  constraint  on  phase  and 
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amplitude  variations  has  been  deduced  by  Kronauer  and  Musa 
in  their  study  of  the  exchange  of  energy  between  oscillations 
in  weakly-nonlinear  conservative  systems. 

In  the  model  of  three  resonant  tanks  with  a  nonlinear 
element  discussed  here,  we  have  derived  a  third  constraint. 
This  constraint  gives  an  additional  relation  among  the  energy 
stored  in  each  oscillator  and,  under  a  special  circumstance, 
it  can  be  reduced  in  the  form  of  the  Manley-Rowe  relation^. 
With  these  three  constraints,  the  behavior  of  three  coupled 
oscillators  can  be  reduced  to  a  problem  of  quadrature. 

We  have  also  extended  the  analysis  to  the  situation  in 
which  there  are  dissipation  elements  in  the  resonant  tanks. 

In  such  a  simplified  model,  the  effect  of  damping  imposes  a 
similar  influence  on  the  energy  exchange  as  that  of  detuning 
of  the  oscillation  frequencies.  If  an  external  source  is 


I 
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provided,  to  make  up  for  the  energy  loss  in  the  system,  a 
steady  state  will  be  reached.  When  the  driving  intensity 
exceeds  a  certain  level,  a  response  with  an  oscillation  fre¬ 
quency  less  than  that  of  the  external  source  will  appear. 


1.  Coupling  of  Three  Matched  Oscillators 

In  a  linear  system  containing  dissipationless  elements 
only,  thpre  are  in  general  some  resonant  modes  and  correspond¬ 
ing  normal  coordinates.  Tranforming  to  these  normal  coordi¬ 
nates,  the  system  may  be  described  by  a  set  of  separated 
second  order  differential  equations* 


2 

d*X 


dt 


+  u)^X.  _ 


l  -  0 


(2-1) 


where  is  the  angular  resonance  frequency  pertaining  to  the 
particular  mode  and  is  the  displacement  in  the  correspond¬ 
ing  normal  coordinates.  Any  disturbances  to  the  system  may 
excite  some  of  these  modes.  The  resultant  response  of  the 
system  is  the  linear  combination  of  these  individual  motions. 

If  the  system  possesses  some  nonlinear  elements  or  the 
disturbance  becomes  sc  strong  that  elements  lose  their  linear 
property,  the  system  will  act  quite  differently*  The  non¬ 
linearity  will  cause  interaction  among  modes,  thereby  negat¬ 
ing  the  linear  superposition  principle.  However,  if  the  non¬ 
linear  effects  are  sufficiently  weak,  we  can  develop  a  per¬ 
turbation  expansion.  The  first  order  expansion  theory  is 
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known  an  the  quasi-linear  approximation. 

In  the  following,  we  will  discuss  a  system  of  three 
resonant  modes  with  a  single  nonlinear  element.  Such  a 
simplified  model  will  illustrate  the  basic  mechanism  of  non¬ 
linear  coupling  without  involving  too  many  mathematical  com¬ 
plications  and  will  provide  a  framework  for  understanding 
the  ohenomena  which  we  shall  pursue  later. 

We  consider  that  the  single  nonlinear  element  is  an 
energy  storage  device  with  the  energy  function  given  bys 


Vn  =  tffXlX2X3  *  (2-2) 

where  e  is  a  nondimens ional  quantity  with  a  magnitude  of 
much  less  than  one,  and  a  is  a  constant  related  to  the  non¬ 
linear  element.  This  might  be  difficult  to  realize  physi¬ 
cally,  but  is,  mathematically,  an  appropriative  device.  The 
total  Lagrange  of  the  three  oscillators  with  such  a  non¬ 
linear  element  is: 


ai*i 

L  =  -P- 


*2 

a2X2 


2  2 


n  n2„2 


+  a3X3.  .  _  .  _  a2“2X_2_  _ 

2  2  2  2 


+  t'  ^ 


(2-3) 


where  a*s  and  w’s  are  parameters  related  to  the  system 
elements. 

The  response  of  such  a  system  is  governed  by  the 
following  coupled  equations: 


\ 
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I 

! 


♦  u)^X1  =  ta1x2x3  , 


— 5^  +  u)tx0  =  ezO-X'A,  . 
dt2  22  23-' 

d2X,  2 

— 5^  +  urx,  =  to,x,x,  . 

at2  ■»  3  3  1  2  ’ 


(2-4 ) 


where  d!s  are  constants  derived  by  the  relation  cf.= 

?m  1  &L^  * 

We  will  at  first  simplify  things  by  assuming  that  the 
o)*s  obey  the  relation, 

+^2  =  ^3  *  (2-5) 

In  Eq*  (2-4)  c*s  are  a  measure  of  the  interaction  among  X^, 
X2,  and  X^*  Without  such  coupling,  X^,  Xg,  and  X-,  will 
oscillate  independently  with  angular  frequencies  to.,  u2, 
and  td~  respectively.  We  refer  to  the  matched  condition  of 
Eq.  (2-5)  as  the  "sum  rule"  for  angular  frequencies  of 
coupled  oscillators. 

Since  eo2,  eCj  are  small  quantities,  we  can 

linearise  Eq.  (2-4)  by  expanding  the  variables  in  terms  of 
6$ 

Xi  "  Xi0  +  «il  +  eZxi2  + . •  i  =  2.  3- 

tf  -  t(l  +  b16r2  +  b2fe3  + . )  , 

tg  =  t(fc+  d^2  +  d263  +  6*'#,,)f 


(2-6) 


1 


where  b's  and  d's  are  constants  of  expansion. 

The  dependent  variable  t  ic  now  expressed  in  two 

scales  t~  and  t  •  We  choose  that  t  is  less  than  t «  by 
is  s  1 

order  of  e ,  thereby  describing  the  response  of  the  system 

in  two  time  domains:  a  fast  time  scale  t^  to  indicate  the 

immediate  response,  and  a  slow  time  scale  t  to  show  aver- 

s 

aging  effects  due  to  the  nonlinear  terms  in  the  equation. 
This  mathematical  approach  has  been  adopted  from  the  two- 
timing  perturbation  method.-^0 

The  differential  operators  then  are  changed  to: 


a  dtf  .  a  dts 

3  at*.  dt  at  at 


ei  + 


d2  _  32  ,dtft2  .  232  .  ■  bf  dtsv  32  fdtSt2 

dt2  9t  2  dt  dt  3t  2^dt 

I  s 

2  2 

=  -K-  +  +  ••••  (2-7) 

at2  atfats  » 

ana  the  X.  *s  are  no w  a  function  of  both  t„  and  t  . 

x  is 

By  substituting  Eqs.  (2-6)  and  (2-7)  into  Eq.  (2-4), 

we  can  collect  together  '.erms  of  the  same  order  as  €.  The 

results  are: 


# » 

For  order: 


+  W1  X10  =  °* 
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v20  2 

T2  +  ^2  X20  =  °' 


iP  +tJ3  x30 


(2-8) 


For  order i 


32X,,  ?  -*2 

~  2~  +  ^1  Xil  =  alX20X30  "  2  3tfats  X10, 


9  Xgi  2 

— 2~  +o)2  x21  *  ®2xiox3o  "  2  atfats  x2o,  (2-9) 


^  X31  2  ^2 

'3t  2  +  ^3  X3i  “  03X1OX2O  “  2  ?tfats  *'30, 


and  so  forth  for  higher  order  terms  of  £. 

The  solutions  of  the  differential  equations  (2-8)  take  the 
formt 


X1C  —  ^1  (cd^tj*  t  £ -^ ) t 

X20  =  R2  cos  ^U'2tf  +  ^2)  ’ 


(2-10) 


m  rji  a  a  «•  f  0  \  ^  ^ 

—  “r  PjJ  * 


where  the  R’s  and  /*s  are  functions  of  t  and  have  to  be 

s 

determined.  In  seeking  solutions  for  Eq.  (2-9),  we  sub¬ 
stitute  the  assumed  solutions  from  Eq.  (2-10)  into  the 
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right  hand  side  of  Eq.  (2-9)  as  source  terms »  After  ex¬ 
panding  terms  X^QX.,0,  Xi0X20,  and  X20**30'  we  ‘fc^a't  they 

contain  the  same  frequency  components  as  the  natural  fre¬ 
quency  of  the  first  order  equations#  This  means  that  the 
X^'s  would  grow  linearily  with  time,  which  would  violate  the 
assumption  that  the  X^  terms  remain  smaller  than  X^Q  terms. 
These  terms,  called  secular  terms,  can  be  avoided  if  we  set 
the  excitation  components  equal  to  zero#  The  process  of 
suppressing  the  secular  terms  gives  the  conditions  for 
determining  the  R*s  and  /f*s.  We  then  obtain  the  following 
relations  from  Eq.  (2-9): 

3R,  CF, 

26)iat ~~T  R2R3  sin  r  =  0, 

3R2  02 

2u)2  dt  2~  R1R3  sin  r  =  °» 


where  r  =  -/2 
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Since  the  R*s  and  /i*s  are  now  functions  of  t  only,  the 

s 

partial  differential  operators  can  be  changed  to  ordinary 
differential  operators.  By  further  combining  the  last  three 
equations  in  Eq.  (2-11),  we  have  in  this  first  order  ap¬ 
proximation  four  coupled  equations  to  describe  the  response 
of  the  system. 

dR-,  ff, 

dt~  -  53  ,K2R3  sin  r  =  °' 

s  1 


dR. 


dt;  -  55,R1R3  sin  r  =  °> 

ra  ,B1R2  sin  r  =  °- 


c 

dts  563 


(2-12) 


and 


dr 

dt 


i*3RlR2 

l*¥T 


®1R2R3  ff2RlR3 

!E¥i  '  '“PT 


cos  r  =  0. 


Eq.  (2-10)  indicates  that  the  response  is  mainly  dominated  by 
the  oscillation  of  the  system  at  its  linear  resonance  frequencies, 
and  Eq,  (2-12)  shows  how  the  amplitudes  and  phases  of  such 
oscillations  are  modified  slowly  in  time  due  to  system’s 
nonlinearities.  If  the  initial  conditions  are  known,  the 
transient  behavior  at  any  instant  can  be  determined  by 
integrating  Eq.  (2-12), 

Prom  the  first  three  equations  cf  Eq.  (2-12),  we  will 

obtain 

n  \2  /-  \  T-,  \2 


“W  ,  „ 

— - -  +  — - -  +  — —  -  c, 


(2-13) 


I 


i 


2-0 


where  E  is  a  constant, 


<"].*?>  (1°2  4>  . 

"l  *  °2 


(2-14) 


where  C  is  also  a  constant. 


A  relation  between  phase  and  amplitudes  can  be  also  be  de¬ 
duced  from  Eq.  (2-12) : 

RlR2H3  cos  r  =  K,  (2-15) 

where  K  is  a  constant.  With  these  constraints  in  Eqs.  (2-13). 
(2-14),  and  (2-15),  the  response  of  the  system  can  be  inte¬ 
grated  out  in  terms  of  slow  time  tg. 

There  are  some  physical  interpretations  associated  with 
the  constraints  of  Eqs.  (2-13),  (2-14),  and  (2-15).  Since 
the  square  of  amplitude  is  a  quantity  for  the  energy,  the 
constant  E  in  Eq.  (2-13)  is  related  to  the  total  energy  ini¬ 
tially  stored  in  the  system.  Eq.  (2-13)  is  just  an  expres¬ 
sion  for  energy  conservation  and  always  exists  for  a  con¬ 
servative  system. 

Equation  (2-14 )  indicates  that  a  certain  relation  has  to 
follow  for  energy  exchange  between  two  modes.  It  depends 
on  the  particular  nonlinearity  we  have  assumed  for  the  system. 
If  the  R's  are  inversely  proportional  to  the  angular  frequen¬ 
cies,  the  w'sand  s' s  are  equal ,  and  the  initial  condition  makes 
the  constant  C  equal  to  zero,  Eq.  (2-14)  will  have  the  form: 


•*1  ^2 


for  Ai=o)iRi 


i=1.2, 3 


(2-16) 
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Manley  and  Rowe^have  derived  this  relation  from  the  property 
of  a  nonlinear  capacitor.  The  constraint  expressed  by  Eq. 
(2-14)  is,  of  course,  implied  to  more  general  cases. 
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According  to  the  analysis  by  Kronauer  and  Musa  ,  the 
amplitude  and  phase  relation  for  a  weakly-nonlinear,  conser¬ 
vative  system  is  determined  by  the  average  value  of  the 
incremental  Lagrange  of  the  system.  As  the  incremental 
Lagrange  in  the  system  discussed  here  is: 

l  =  eax^Xy  (2-17) 

the  average  value  of  i  can  be  found  by  using  the  results  of 
Eq.  (2-10):  ,  T 

l  =  -i-  ««W3  dv  (2-18i 

Jo 

—  (rffR^RgR-j  . 

Hence  the  constraint  of  Eq.  (2-15)  is  due  to  the  incremental 

~0 

Lagrange  of  the  system  with  K  =  — ,  a  general  result  of 
first  order  approximation. 

The  transient  behavior  of  such  coupled  oscillation  can 
be  graphically  illustrated  by  a  topological  approach  as  in 
Fig.  2-1.  The  diagram  uses  three  axes  to  express  the  mag¬ 
nitude  of  fcijR-L,  cOgRg,  an(*  .  Any  point  in  the  space 
represents  a  instantaneous  state  of  the  system’s  response, 

A  surface  of  constant  E  is  an  ellipsoid  on  this  space.  0^, 
Og,  and  are  singular  points  and  they  are  stationary 
states.  However,  0^  is  an  unstable  stationary  state.  Any 
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disturbance  around  0^  will  cause  a  response  with  a  motion 
locus  circling  around  the  ellipsoid.  This  indicates  that 
the  energy  associated  with  the  high  frequency  mode  is  unstable 
and  can  easily  be  diverted  to  the  energy  in  lower  frequency 
modes.  Systems  exhibiting  such  features  can  be  excited 
into  subharmonic  oscillation.  Wt.  shall  further  explore  such 
phenomenon  later. 

2,  Effect  of  Detuning  on  Coupled  Oscillation 

V/e  now  consider  the  case  when  the  frequencies  of  the 
three  oscillators  are  not  exactly  matched.  The  deviation 
from  Eq.  (2-5)  can  be  accounted  for  by  solving  a  detuning 
problem.  We  can  expand  the  unmatched  relation  in  terms  of 
the  matched  condition  as  follows* 
as  o)L  +  /  tdv 

we  define  + 

W2  ~  °°Z  +e62’  (2-19) 

0O3  =  * 

and  let  cOJ  +  =6 Oj, 

where  6’s  are  detuning  parameters  from  their  corresponding 
modes. 

The  original  coupling  equations  can  be  rewritten  as* 
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~  +  to*  ?XJL  =feCX2X3  +  2fe61u)|X1 


+  cO’?X~  =£<JX,  X-i  + 


*2^2X2 


(2-20) 


2  +  oO'?X3  =^X1X2  +  2£63<X>*X-  , 


where  we  let  all  0‘s  be  the  same  to  simplify  the  analysis. 
Using  again  the  two-timing  perturbation  method,  we  obtain 
the  first  order  solution: 


X10  =  R1  cos  ^lJitf  +  6*1) 

X20  =  R2  co<5  ^^f  +  ^2^ 
and  X^q  =  R0  cos  (o^t^  + 

with  the  auxiliary  relations: 
dR,  aR?R^ 

at  r  °  • 

s  l 


d R~  ffR.rt- 

dt  sin  r  =  0  ? 

s 


dR-j  ^*1^2 

dtl  +  "^T~  3ln  r  "  c  • 


(2-21) 


(2-22) 


Hi"  R->8?  RoR-5  R.R:, 

dtT  +  (  Tuojp^  °  cos  r  +  6  =  0  ’ 


•where  6  =  6^  +  62+6^. 

Comparing  witr.  Eos,  (2-10)  and  (2-12),  we  find  that  the 
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total  detuning  factor  6  appears  in  the  phase  equation  only. 
We  shall  analyze  its  effect  through  the  motion  locus  in  the 
phase  diagrams. 

Changing  R's  into  new  variables,  we  get: 

A1  =toiRi'  A2  =  to2R2’  A3  =  ^3^3*  (2-23) 

The  corresponding  constraints  become: 

A1  +  A2  +  A3  =  E  -  e?.  (2-24) 

A2  A2 

A-.  A-j  ^ 

*77  -  ;tt  —  0  —  2.  a  if.  (2-25) 

CO*J  W£ 

A.,A2A3  cos  r  +  26co^u)'A^  =  K,  (2-26) 

where  we  have  defined  new  constants  e  and  ft. 

Since  these  three  constraints  are  directly  derived  from 
Eq.  (2-22),  the  original  four  variables  A^,  A2,  A y  and  r 
for  describing  the  state  of  the  system  in  four-dimension 
space  can  be  reduced  to  quadrature  and  the  time  variable  car. 
be  found  by  integrating  out  such  a  relation. 

The  ellipsoid  diagram  of  Fig.  2-1  is  still  appropriate 
to  represent  the  response  of  the  system,  but,  because  of 
detuning,  it  is  no  longer  possible  to  deduce  the  phase,  r, 
unambiguously  from  that  diagram.  That  is,  detuning  makes 
possible  monotonic  phase  change  (  while  with  6=0  phase 
oscillates  periodically  ).  To  observe  the  phase,  we  make 
a  transformation  based  on  the  first  two  constraints,  Eqs. 
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(2-24)  and  (2-25),  introducing  a  new  variable  OC  defined  by: 


A 


1 


0>?  p 

e  (  sin  06  + 


K 


)* 


f 


and 


(op  o  i 

Ao  =  e  (  -4  sin  X  -  K  )a, 

kj  «  e  cos  x  • 


(2-2?) 


The  variable  X  is  a  measure  of  the  extent  of  energy  exchange 
from  A^  to  the  and  A2  set.  The  third  constraint  now 

becomes: 


3/^1 ^  _;„4  ~  ,.2  X  _  .  ..  .. 

e  (  ,-j  sin  X  —  -y  ?( sin  X  —  K.  )  cos  X  cos  r 

+2e26u)J4£  cos2X=  K  .  (2-28) 

The  properties  described  by  Eq,  (2-28)  can  best  be 
understood  by  considering  four  cases.  First,  we  observe 
that  the  parameter  &  represents  a  fixed  energy  imbalance 
between  oscillators  1  and  2  and  is  unaltered  by  an  exchange 
of  energy  between  this  pair  and  oscillator  3«  The  parameter 
6  of  course  represents  detuning.  Fig.  2-2  shows  the  solution 
trajectories  (  plots  of  Eq.  (2-28)  for  various  fixed  values 
of  <K  )for  the  four  interesting  cases. 

Case  I.  (Fig.  2-2a).  With  no  detuning,  the  trajecto¬ 
ries  remain  in  a  cell  of  width  it.  Since  the 
imbalance  is  zero,  the  full  rang;  cfXis 
accessible . 


Case  II.  (Pig.  2~2b),  With  detuning,  same  trajectories 
can  exhibit  steady  phase  change*  Since  the 
imbalance  is  zero,  the  full  range  of  %  is 
still  accessible,  although  there  is  no  single 
trajectory  which  makes  the  swing  from X=o  to 
X--n/2»  Thus  we  can  say  that  the  presence  of 
detuning  acts  to  reduce  the  extent  of  energy 
exchange  between  oscillators  3  and  the  1,  2 
set. 

Case  III. (Fig.  2-2c).  Vfith  imbalance,  the  accessible 
range  of  %  is  reduced.  The  absence  of  de¬ 
tuning  permits  trajectories  with  maximum  ex¬ 
change  (subject  to  the  imbalance  reduction). 

As  in  Case  I,  the  trajectories  are  confined 
to  a  cell  of  width  tt. 

Case  IV.  (Fig,  2-2d).  With  both  detuning  and  imbal¬ 
ance,  two  effects  are  seen  to  reduce  the 
extent  of  energy  exchange  throughout  the 

.  Steady  phase  change  is  also  common. 


field 
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3«  Forced  Excitation  for  Subharmonic  Oscillation 

For  systems  containing  damping  elements,  a  transient 
disturbance  will  not  excite  and  sustain  subharmonic  oscil¬ 
lation,  because  the  signal's  limited  energy  will  be  quickly 
dissipated.  In  general,  in  order  to  sustain  a  steady  re¬ 
sponse,  a  constant  energy  source  should  be  provided  to  com¬ 
pensate  for  the  energy  loss  in  the  system. 

From  our  study  of  the  transient  behavior  of  a  conser¬ 
vative  system,  we  have  learned  that  the  energy  in  the  high 
frequency  mode  can  be  easily  converted  into  the  energy  in 
lower  frequency  modes  by  some  disturbance.  We  shall  now 
investigate  the  criteria  for  sustaining  subharnionic  oscil¬ 
lation  when  a  constant  energy  source  is  applied  to  the  system. 
We  use  the  same  model  for  three  coupled  oscillators  but 
include  a  damping  element  in  each  oscillator.  The  original 
three  oscillation  frequencies,  u)^,  cu-g*  are  not  exactly 
matched,  that  is  The  external  source  with 

angular  frequency  ujq  has  excited  two  su'oharmonics  of  angular 
frequencies andu)£  with  the  following  relatiom 
c«Jq  =  t»->  *  +  , 

uli  =c^i  +  e63  f 

u>2  =C02  +  £62  _  (2-29) 

where  6*s  are  detuning  factors  for  each  oscillation  mode  and 


2-19 


6  is  a  small  constant  of  magnitude  less  than  one. 

The  mathematical  formulation  of  this  problem  is  ex¬ 
pressed  in  Eq.  (2-30): 

d2X-  dXq  , 

— ^  +  <:»_ — 2-  +  totrX-,  =  6cr.,X,X0  +  2e;ai-X-J 
dt2  (3dt  0  3  3  1  2  0  3  3 


0"  3  3 


~eF  sin  coQt  , 

d2X,  dX,  „ 

_i  +  61^-i  ♦  .  £aiX2X3  +  2«UJilXl  . 


and  A? 


~"u  — ^  +  en0 — -  +o;«2XC)  =  e<30X,X^  +  26UJI6  X, 

dx  *dt  *  *  c  1  '  t  t  , 

(2-30) 

where  we  have  considered  that  the  dissipation  factors,  »j*s, 
coupling  coefficients,  a's,  and  the  forcing  term  due  to  the 
external  source,  F,  are  small  order  terms. 

We  again  use  the  two-timing  perturbation  method,  expanding 
the  dependent  variable  t  in  a  fast  scale,  tf,  and  a  slow  one, 
t_.  After  summing  up  functions  of  same  order,  we  have; 

O 

for  e  order: 

772  X30  +U)0X30  =  0  * 

dtj. 


2  X10  +aJl  X10  *  0  * 


(2-31) 


2  X20  +UJ2  X20  =  0  ? 


and 
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for  6  orderi 


at2X31  +  toQ2x31  =  '*23-tf3tsX30  +  t?3X10X20  *  2<°063X30 


‘^3H^X30  ‘  F  sin  "o*  • 


2X11  +a){2xil  =  "2at.3t  X10  +  alX20X30  *  2o;i6lX10 

I  s 

-h^ho  , 

2  2 

and  “jX21  +W22x21  =  X20  +  tf2X10X30  +  2tO262X20 

jl  S 


(2-32) 


The  solutions  of  Eq„  (2-31)  aret 
X^q  =  R^  cos  ^ 

xio  =  Ri  cos  + 

and  X2Q  =  R2  cos  \W|tf  +  /f2) 


(2-33) 


After  eliminating  the  secular  terms  in  Eq.  (2-32),  we  find 
the  supplemental  equations  for  solutions  in  Eq,  ( 2**33 ) • 

3R-, 

2uJbat  +  ^3^R3  +  2^rir2  t:in  r  “  F  cos  ^3  =  0  * 

S 

^iat-  *  ^i^i  "  2^R2R3  sin  r  -  0  * 


aR?  c? 

2u?2at~  +  72^2R2  "  2“R1R3  sm  r  =  0  , 
s 
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a/fo  ffo 

2woR3at  +  2  rir2  cos  r  +  2a)o63R3  +  F  sin  =  0  » 

S 

a/,  cr, 

2cO'Ri~t  +  2iR2R3  cos  r  +  =  0  »  (2-34) 

S 

9/^2  &  o 

and  2t02R2aF  +  2rR3  C0S  r  +  2a;262R2  =  °* 
s 

For  transient  behavior,  we  let  F  =  0.  The  total  energy  of 


the  system  to  zeroth  order  is: 

(0«R,)2  (cO*Rp)2  HR  ): 
E  =  — —  +  — ~ —  +  — —r2— 
°1  °2  °3 

By  the  relations  in  Eq.  (2-3^),  we  have: 


.  (2-35) 


is  “  -  57<“'iri)2  -  !f<“’2R2)2  -  ^Hr3)2  • 


(2-36) 


It  is  convenient  to  let  7l=  ?2=  ^y=^it  ^q.  (2-36)  becomes: 


=  -  or  E  =  EQe' 


(2-37) 


This  indicates  that  the  total  energy  will  eventually  be  dis¬ 
sipated.  The  ellipsoid  of  Fig.  2-1  will  shrink  to  zero  in 


one  xj-iujLo  ujl  j.Ong  UJLSIIC  • 


The  forced  response  can  also  be  expressed  in  terms  of 
the  system’s  energy  by  keeping  the  F  term  in  Eq.  (2-36)*  The 


result  is: 


-  ■  V 


Fu^R, 


cos  /- 


(2-38) 
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There  is  a  competition  between  the  energy  dissipated  in  the 
system  and  the  energy  supplied  from  external  source,  and  a  sus¬ 
tained  driven  oscillation  will  correspond  to  a  balance  of 
this  effect. 

Another  interesting  feature  of  this  nonconservative 
system  is  that  the  previous  constraint, 


C  = 


l“l 


er. 


(2-39) 


W1  2 

is  no  longer  a  constant.  By  the  relations  in  Eq.  (2-3*0. 
we  can  determine  its  variation  ast 


dC  7ia>iRl  72u)2R2 

"  ■;  ■  —  «•  ■■■■'  —  —  ■—  ‘f  ■—  " 

dt„  CJo 


(2-40) 


For^1=72=^=,I,  ^  can  further  simplified  to 

-P  -tft 

^  =  -  rjC  or  G  =  C0e  s  .  (2-41) 

a  s  *  u 

For  the  steady  state  (t..— ►»),  C  will  become  zero  in  this 
case  of  balanced  dissipation.  Hence  Eq,  (2-4l)  gives  the 
exact  condition  for  the  Jfianley-Rowe  relation^. 

The  response  of  each  frequency  component  after  the 
system  reaches  a  steady  state  car,  be  worked  out  by  letting 
all  the  time  variable  terms  in  Eq*  (2—34)  go  to  zero j  A 
condition  for  detuning  factor  8's  is  obtained  ast 


a2  *7  2 


(2-42) 


or  we  can  define  a  new  constant  A  as  the  ratio  of  6  to  ^ 


s 
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in  the  following  way* 

6,  6? 

A=r^  =  ~  (2-43) 

’ll  1z 

Eq.  (2-43)  shows  that  the  detuning  is  proportional  to  the 
dissipation  factor  and  it  shifts  in  same  direction  according 
to  whether  the  driving  signal  frequency  is  larger  or  less 
than  the  sum  of  two  original  subharmonic  frequencies. 

The  amplitude  of  each  frequency  component  is  given  by: 


2  „  ““Hfe 


R,  =  „  „ 

1 


2  3 


[' 


j(^3  -  4Afi3)2  +(1  4  4A2) 

-  ( 3  ^  “  ^3  “^3^1  * 


?2°1Q2 


|(^3  -  IJA63)2  +(1  -f  4A2) 


2  _  M^o7i 

2  o,o 


13 


and  R3  = 


2  H^STlfe 


12 

«ltf2 


(  1+4^) 


( 2—44 ) 

Since  all  R*s  are  physical  quantities,  they  have  to  be 
real  values.  The  minimum  intensity  of  P  for  R^  and  Rg  to 
appear,  or,  in  other  words,  the  threshold  for  generating  the 
subharmonic  components,  is  determined  by» 

02  _  /K*V  2  011^2  yjZ  +  4 ^)(X  +  4A2)  .  (2-45) 


th 


<*la2 
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Below  this  threshold,  the  only  significant  response  is  R^. 
which  is  described  by  i 


P 

cdQ(^  +  462)? 


(2-46) 


We  see  that  exhibits  a  very  interesting  behavior*  it  in¬ 
creases  linearly  with  the  external  driving  force  and  then 
stays  at  a  constant  value  after  the  threshold  for  subharmorJc 
generation  is  reached.  Pig.  2-3a  is  a  typical  graphical  dis¬ 
play  of  the  change  in  R^,  Rg,  and  R^»  with  respect  to  the 
increase  in  the  intensity  of  P.  Pig.  2-3b  shows  the  experi- 
mental  data  collected  by  Y.  Tsuzuki  and  M.  Kakuishi  in 
their  observation  of  the  excitation  of  contour  modes  of  vi- 
braxion  in  AT-cut  quartz.  The  responses  of  the  subharmonics 
have  been  squared  in  that  diagram,  so  they  appear  to  be 
linear  past  the  threshold.  They  fit  our  mathematical  model 
and  analysis  very  well. 
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FIG.  2-3  RESPONSE  OF  FORCED  EXCITATION 
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Chapter  III 

SUBHARMONIC  GENERATION  IN  DISTRIBUTED  SYSTEMS 

The  effect  of  inxense  acoustic  waves  has  been  studied 

by  various  authors,  for  example t  Ghiron*^,  Hunt-^,  and 
'34 

Andreev-^  ,  They  indicated  that  the  nonlinear  properties  of 
the  medium  as  well  as  convection  of  the  disturDance  usually 
cause  wave  distortion  during  propagation  arid  the  generation 
of  superharmonics  of  the  original  signal.  Experimental 
observations-^  have  confirmed  these  analytic  predictions. 
However  such  a  theoretical  approach  has  failed  to  supply 
an  explanation  for  the  phenomenon  of  subharmonics.  We  shall 
provide  a  detailed  study  of  this  aspect  of  the  problem. 

In  order  to  simplify  the  mathematical  operations,  the 
entire  analysis  throughout  this  chapter  is  based  on  a  one¬ 
dimensional  model.  Euler ian  coordinates  are  employed  in  all 
derivations  to  facilitate  comparison  of  the  theoretical  re¬ 
sults  with  the  laboratory  data  from  the  experimental  work 
that  is  described  in  the  next  chapter. 

We  review  the  formulation  of  the  acoustic  wave  equation 
from  the  basic  conservation  laws  in  the  beginning  of  this 
chapter*  As  the  primary  interest  of  this  study  is  to  look 
for  the  essential  contribution  to  the  subharmonic  generation 
for  an  intense  acoustic  wave,  the  order  of  magnitude  of  all 
pertinent  parameters  in  the  system  considered  is  examined. 
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The  two-variable  method  is  adopted  to  fir1  an  approximate 
solution  for  the  nonlinear  wave  equation.  We  obtain  the 
energy  exchange  relations  among  three  signals  when  their 
frequencies  follow  a  sum  rule. 

This  analysis  is  then  extended  to  a  multirasonator, 
and  the  threshold  of  exciting  a  subharmonic  pair  is  deduced. 
We  also  discuss  the  effects  of  detuning  and  energy  dissi¬ 
pation  on  subharmonic  generation* 


1.  Formulation  of  Finite  Amplitude  Acoustic  Wave  Motion 

The  basic  laws  describing  acoustic  wave  motion  can  be 

derived  from  the  hydrodynamic  equations.  In  one-dimensional 

* 

Eulerian  form,  they  are  i 

Conservation  of  mass 


3P*  .  afu’  _  n 

at*  ax’  "  u 

Conservation  of  momentum 

3f*u 


(3-1 ) 


_L  +  u'a»£l  +  ap*  -  o  (3-2) 

at*  +  uax’  +  ax’  »  _  .2  “ 

O  A 

where  p’  and  u*  are  the  density  and  the  particle  velocity 
respectively  and  m!  ip  the  coefficient  of  viscosity  which 


n  f\r> ci  o+a  ryP 


coefficient  of  shear  .viscosity  \j'  arid  ihe 


coefficient  of  dilatation  viscosity  X'  ,  as  defined  by  the 


*  The  primes  are  used  to  denote  dimensional  quantities.  We 
will  shortly  reduce  them  to  nondimensional  variables  and 
the  primes  will  be  dropped. 


relation* 

=  2V*  +  X‘.  (3-3) 

The  losses  due  to  viscosity  are  small,  and  so  are  the 
losses  due  to  heat  conduction.  Under  these  conditions,  the 
process  can  be  assumed  to  be  adiabatic  as  far  as  wave  speed 
is  concerned  since  the  losses  affect  the  speed  only  quad- 
ratically.  The  state  variable  p*,  the  pressure,  can  be  ex¬ 
pressed  by  expanding  around  the  equilibrium  state  as  follows* 
P'  =  p'(e*). 

=  P°  +  +  * . 

=  v'0  +  c  *oUee)  +  |  r’^e*)2  +••••••••*.  (3-*0 

c^  is  so  called  the  small  signal  sound  speed,  p  *  is  a  non¬ 
linear  parameter  and  Af*  =  £,-€q  is  the  deviation  of  the  den¬ 
sity  from  the  equilibrium  state.  The  subscript  q  denotes 
that  the  expansion  is  assumed  adiabatic. 

Some  authors-^' in  discussing  the  nonlinear  be¬ 
havior  of  a  fluid,  write  Eq ,  (3-^)  in  the  form 


(3-6) 


'*-h 


Tn  the  case  of  an  ideal  gas,  these  coefficients  can  be  re¬ 
lated  to  the  ratio  of  specific  heat  Y  * 

A  =  YF^,  B  =  Y (Y-l )t>q»  and  I  =  Y-l.  (3**?) 

'By  using  the  approximate  relation  of  Eq.  (3-^)»  Eqs. 
(3-1)  and  (3-2)  then  can  be  written  in  terms  of  <o*  and  u' as 
follows! 


ae'u1 


3  xc 


0, 


+  u 1 


ag’u*  +  c.2  a(Aey). 


ax' 


ii 

a,-2 


(3-8) 


In  order  to  make  an  approximation  with  regard  to  the 
order  of  magnitude  of  various  terms  in  Eq.  (3-8),  we  put 
Eq.  (3-8)  in  nondiraensional  form  by  adopting  the  following 
scaling  parameters! 

Independent  variables! 


Length  •  •  . . .  X  cm 

Time  »••**••••••  T  sec 


Dependent  variables!  (3-9) 

Partical  velocity  ••••••••  u  cm/sec 

Density  ••••••••••••••••••  gm/ctP 


Then  Eq.  (3-8)  has  the  formi 
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ae  ut  a£u  =  o 

St  X  Sx 


agii  +  Ol  +  i 

^  O  U  "?N  V  *  O  -VV  “  O 


r'  t2p 


x  at 


^2  bx 
2  2 

JU?UT  S  u 

2  2 
X^e0  ax* 


=  o. 


We  denote  the  nondimensional  constants  byt 


3(Ag)2 

S  x 


(3-10) 


"  B  “X  • 

n-rif£t 

I  "  9 


(3-U) 


= 


c  •  2m2 

c  o1 


(1  —  . 

^  e0XU  * 

Eq.  (3-10)  then  becomes* 


+  M  a|u_  =  0, 


e  *  A  fr  ♦  «2  iff- ♦  4r*gfi-  - 


=  0. 


(3-12) 


For  the  case  of  a  propagating  harmonic  wave  where 
both  losses  and  nonlinear  effects  are  small,  the  natural 
scaling  parameters  to  use  are* 


(3-13) 


whereto^,  and  are  the  dimensional  angular  frequency  and 


Then 
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9 

(Or 

wave  number  of  the  wave, respectively,  and  r- r~  =  cl. 

Kr  v 

UKr  U 

M  =  — r  =  ~  is  seen  to  be  the  particle  velocity  Mach 
<°r  c0 


number 


and  c  =  cl(— )  =1,  The  problem  of  interest  here  is 
u  oJU 


when  M  is  small,  but  not  negligibly  so.  Therefore  we  intro¬ 
duce  the  symbol 


6  =  M  = 


U 


(3-1*0 


and  Eq.  (3-12)  becomes 


§£  +  ^  aeu  _  Q 
dt  *  6  9x  "  °* 


+  aa  ♦  3<ae)-  +  1  r2  -Stegl!  .  e2^-3^-  *  0. 
ax  ax  ax  2 1  ax  e  ^  ax2 


Ur 

e  =  m  =  « 1. 

co 


(3-15) 


2,  Wave  Equations  in  Terms  of  Slow  and  Fast  Variables 

To  get  the  explicit  expression  for  e  and  u  as  a  func¬ 
tion  of  x  and  t  by  solving  Eq.  (3-15)  is  a  rather  difficult 
mathematical  task.  An  approximate  method  can  be  developed 
for  the  case  of  a  small  Kach  number.  First,  consider  the 
linearized  invi-cid  form  of  Eq.  (3-15).  (Remember,  density 
is  nondimensionalized  so  1  at  equilibrium). 
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at  ax  u‘ 


au 


(3-16) 

at  ax  Uf 

which  further  reduces  to  the  form  of  a  wave  equation. 


o 

a  u 


=  0. 


(3-17) 


ax"  at 

The  corresponding  solution  of  Eq,  (3-1?)  has  the  form 
u  =  f(x-t)  *  g(x+t),  (3-18) 


where  f  represents  tne  wave  travelling  to  the  right,  a  for¬ 
ward  wave,  and  g  represents  the  wave  travelling  to  the  left, 
a  backward  wave. 

For  a  simple  sinusoidal  excitation  with  an  angular 
frequency  a), 

u  =  A  sin  (cot-kx+/f)  +  Bsin  (uot+kx*#),  (3-19) 

where  A,  B,  /f,  and  $  are  constants  which  depend  on  the 
initial  condition  of  excitation. 

Eq,  (3-19)  describes  a  wave  propagating  with  a  con¬ 
stant  amplitude.  For  finite  amplitude  wave  motion,  if  e 
is  much  less  than  one,  the  physical  picture  of  Eq.  (3-15) 
should  not  be  very  different  from  Eq,  (3-16),  The  convec¬ 
tion,  nonlinear  parameter  and  viscosity  may  just  modify 
slightly  the  propagating  wave  form.  Since  their  effect  is 
small,  it  may  be  considered  that  the  original  wave  form  has 
a  slow  variation  in  the  amplitude  and  phase,  with  respect  to 
time  and  space,  as  it  propagates  through  the  medium.  With 


3-8 


bles  are i 


3-9 


at  <9t.p 

a  _  a 


5  +  6  5 


at. 


ax  ax 

J2. 


+  £ 


axs  ' 


(3-22) 


and 


2  2 
ax  ax^ 


a  +  26  a 


axfsxs 


+  & 


2  a£ 


ax 


2  • 


Substituting  Eqs.  (3-21)  and  (3-22)  into  Ea.  (3-15) 
we  have i 

+e^~)(1  +*ei  +  ^^2  +  +6fe) 

is  is 


*  (l  +  ee1  +  e  e2  + 


)(ul  +  tu2  + 


+  e<?i +  *2e2 + 


•  •  •  )  =  0, 

)  (u^  +  ^U2  . . ) 


+  e2(ux  +eu2  . . )(—;  +  *  Jr>(  1  +  «<?i  +  ♦••••> 


(u^  +  eu2  + 


> +  <ife+e#:>(ee  1 +  £2<?2+-— 1 


+  +  &2ea  + 


)■ 


-  <M(^2  +  +  e2^2)(ul  +  £u2  . . >  =  °' 


ax. 


ax. 


(3-23) 


In  Eq.  (3-23).  the  order  &1  yields 


“at-^i  +  ax-  ui  “  °* 


atf  ui  +  sxf^i  "  °’ 


(3-24) 


! 
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v/hich  can  be  reduced  to  a  first  order  wave  equation  by 
eliminating 


u„  - 


ax 


2 


u.,  =  0. 


at 


2  “1 


(3-25) 


For  order  e  »  we  have  the  following  equations i 

atj^2  +  atT^l  +  u2  +  ax“  ^lul  +  3X^  u2  =  °' 


0  3  ^  3  a  a  o 

atf  u2  stf  e  iui  atg  ui  uiaxf  ui  s>x~  ^  2 

+  +  -  u,  r  o. 

5  I  C*X  jp» 


(3-26) 


The  corresponding  second  order  wave  equation  is: 


ax. 


7^  U2  "  at2  u2  ~  _2  elul  “ 


at; 


e1u1  +  ~- 
11 


luiJb  ui) 


o  3  o 

.  r  a^  „2  >11  .  s 

2  atfaxf  e  l  *  2  _2^.  ui  *  at^t„  ui 


ax^at^ 


'f ^  s 


"  3xf9xs  ui  +  at^ax^  v  i  ax^atg  ^  i' 

(3-2 7) 

The  first  order  solution  can  be  found  from  Eq.  (3-25)*  It  is 
then  used  to  generate  source  terms  in  the  right  hand  side  of 
Eq,  (3-2?),  The  higher  order  solutions  can  be  found  in  a 
similar  way.  However,  this  only  gives  the  general  form 


e,. 


3-11 


of  the  solution.  For  a  particular  problem,  the  initial  and 
the  boundary  conditions  are  required  to  select  the  meaningful 
solution.  In  the  following  sections,  pa^  oicular  cases  of 
subharmonic  generation  will  be  studied  in  some  deteil. 


3.  Interaction  of  Waves 

For  certain  initial  disturbances,  the  response,  in 
terms  of  the  particle  velocity  of  the  medium,  according  to 
Eq.  (3-25),  has  the  general  form: 


u  =  2  T  A  cos£„  +  B  cos  £  1 
n  Ln  ^n  n  ^nj  , 


(3-28) 


where 


An  =  VvV' 
Bn  *  Bn**s,ts^’ 


=  “nV  -  Vf  +  4>(VV* 

=  <Vf  +  krxf  +  *n<VV‘ 


(3-29) 


Suppose  the  disturbance  at  x  =  0  is  mainly  dominated 
by  three  distinct  signals  v/ith  angular  frequencies  oo^,, 


and  cOt  related  by: 

j 


OJ^  +  CJg  =  » 

Then  Eq.  (3-28)  can  be  rewritten  as 

3  r 


(3-30) 


*n  =  n=l  [An  cos§ n  +  Bn  cos^n]  • 


with  a  corresponding  density: 
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=  n=l  [An  cos  h  n  "  Bn  cos$n]  •  (3'32) 

For  nondispersive  media,  the  wave  number  k’s  will  follow  the 
same  sum  rule  as  that  of  angular  frequency  from  Eq.  (3-30) » 

kx  +  k2  =  k3  .  (3-33) 

An*  Bn’  ^r\*  and  are  on3-y  £iven  f°r  ‘t^e  initial  condition 

x  =  0*  "t  =  0.  To  determine  how  they  change  with  respect 
s  s 

to  time  and  space  (in  slow  scale)  depends  on  further  informa¬ 
tion  from  the  second  order  wave  equation  in  Eq.  (3-2?). 

When  we  substitute  Eqs.  (3-31)  and  (3-32)  into  the 
second  order  wave  equation,  Eq,  (3-2?),  to  solve  for  the 
second  order  term  of  u,  v/e  find  that  the  source  terms  on  the 
right  hand  side  of  that  equation  contain  the  same  frequencies 
as  the  natural  frequencies  on  the  left  hand  side.  This  means 
that  the  solution  will  have  terms  increasing  linearly  with 
time  and  space.  Such  solutions  will  ultimately  violate  the 
original  requirement  of  6u2«u^  in  the  perturbation  expansion. 
Those  terms  which  cause  such  limitations  on  perturbation  ex¬ 
tension  are  called  the  secular  terms.  Since  the  amplitude 
and  phase  are  assumed  to  be  functions  of  the  slow  variables 
in  the  case  discussed  here,  the  problem  can  be  avoided  by 
choosing  the  slowly  varying  amplitude  and  phases  in  order  to 
make  the  secular  terms  vanish.  This  choice  leads  to  the  re¬ 
lation  showing  how  amolitude  and  phase  change  with  time  and 
space  (in  slow  scale)  after  the  initial  disturbance. 
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After  some  tedious  algebraic  operations,  we  achieve 
the  following  results* 

3A-  9A1  2 

aFT  +  W  +  ^W1A1  "  tftUlA2A3  sin  r  =  °* 

S  S  ‘ 

3 Ap  ^ Ap  p 

_  +  _  +  >w2A2  -  tfw2A1A3  sin  r  =  0, 
s  s 


3A-  9A-  P 

+  >1W3A3  +  o  w3A1A2  sin  r  =  0, 

Ai  it:  +  Ai  jr  +  0<uia2a3  =°s  r  =  °- 

s  s 

9/p  3^2 

A2  3rt^  +  A2  +  ato2AlA3  C0S  r  =  °» 

A3  atf  +  A3  axf  +  a c°3A1A2  003  r  =  °- 


(3-34) 


9B-  SB..  „ 

5T  "  53T  +  -  ow1b2b3  sin  s  =  0, 

s  s 

9B?  3Bp  P 

at:  -  5x1  +  Mwj2b2  -  0aJ2BlB3  sin  3  =  °- 

s  s 


aBo  3Bp  2 

at  "  +  ^uu^3B3  +  ^WpB^Bg  sin  s  =  0, 

s  s 

Bi  5C  ‘  Bi  5^  +  owl32b3  '=os  s  *  °- 

s  s 

as? 

B0  -^3^  -  B0  - —  +  cos  s  =  0, 

2  atg  2  $xs  213 

3*3  .  a* 

B3  5t  '  S3  +  oWjBjBj  oos  s  =  0, 
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where 

*  *  5  (1  +D  , 

r  «  /rf3  .  ^  4  ,  (3-35) 

s  -  §3  -  -  §2  • 

In  Eq.  (3-34),  the  first  six  equations  are  independent  of 
the  last  six  equations.  This  indicates  that  there  is  no 
coupling  between  An,  /n,  and  B  ,  <5  .  Since,  from  Eq,  (3-29). 

Aft  and  are  amplitudes  and  phases  respectively  of  the  for¬ 
ward  travelling  waves  and  Bn  and  are  amplitudes  and 
phases  respectively  of  the  backward  travelling  waves,  the 
results  expressed  in  Eq.  (3-3*0  show  that  these  two  types  of 
wave  propagate  independently  without  interaction, 

V/e  now  fui  l her  examine  the  charactc  r  of  the  forward 
wave,  to  understand  the  mechanics  of  wave  interaction  (the 
conclusion  will  also  apply  to  the  case  of  backward  wave). 

The  first  six  equations  in  (3-34)  have  a  single  set 
of  characteristics  in  slow  space  and  time  with  slope  unity. 
Therefore  if  we  introduce  the  variable 

r=ts-xs.  0-36) 

these  equations  may  be  written  as» 

dA.j  2 

rpp  +jvU)^a^  -  <3rWlA2A3  s*n  r  =  °» 
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dA, 


-~r  +  }kto  2^2  “  <Jw2A1A3  sin  r  “  °* 


dA 


u  T  2 

+  ju.w^A^  +  ffco^A^Ag  sin  r  =  0, 


3  3 


d^i 


A1  dcf  +  a<JJiA2AJ  cos  r  “ 
d^9 

A2  dc^  +  tfa>2AlA2  cos  r  “ 
d/. 


(3-37) 


A~  +  doo^A.A.,  cos  r  =  0. 

J  Ql  j  j.  £ 

Although  amplitudes  and  phases  change  slowly#  it  is 
interesting  to  see  that  the  slow  behavior  propagates  at  the 
same  speed  as  the  waves  themselves.  Eq.  (3-37)  can  be 
further  reduced  to 

dA.  2 

-  Cix^AgA^  sin  r  =  0, 

dA?  tp 

dr~  +  >“2*2  -  «w2aia3  sin  r  =  °' 

dAo  a 

'rzr-  +  jx<dtk~  +  auu^A^Ag  sin  r  =  0, 


(3-38) 


dz  ’  ^33 

dr  +  (U3AiA2  _  ColA2A3  _  *°2A1A3 
d£  A^  A^  A^ 


cf  cos  r  =  0. 


For  a  nonviscous  medium,  these  coupling  equations  can  be  even 
simplified  toi 


dA. 

-  aco^ AgA^  sin  r  =  0, 

d.A2 

jjj—  -  d c*j2A1A3  sin  r  “  °» 


3-lfe 


where 

«  *  |  (i  +  r)  . 

r  =  ff3  -  ^  -  /f2  ,  (3-35) 

s  ~  $9  . 

3  ».  2 

in  Eq.  (3-3**  )#  the  first  six  equations  are  independent  of 
the  last  six  equations#  This  indicates  that  there  is  no 
coupling  between  Afi,  /n,  and  En,  ®n#  Since,  from  Eq.  (3-29), 

A.fi  and  6n  are  amplitudes  and  phases  respectively  of  the  for¬ 
ward  travelling  waves  and  and  are  amplitudes  and 

n  n 

phases  respectively  of  the  backward  travelling  waves,  the 

results  expressed  in  Eq.  (3-3^)  show  that  these  two  types  of 

wave  propagate  independently  without  interaction r 

We  now  further  examine  the  character  of  the  forward 

wave  to  understand  the  mechanics  of  wave  interaction  ( the 

conclusion  will  also  apply  to  the  case  of  backward  wave}# 

The  first  six  equations  in  (3-3*0  have  a  single  set 

of  characteristics  in  slow  space  and  time  with  slope  unity. 

Therefore  if  we  introduce  the  variable 

f  =  t  —  x  , 

s  s 

these  equations  may  be  written  asi 
dA*  9 

+/v^^A^  -  Oco^ApA^  sin  r  =  0, 


(3-36) 


3-3  5 
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A, A  A„  cos  r  =  K. 

i.  t  ) 

The  mathematical  sisnil  ariliiy  between  wave  interaction  and 
couoled  oscillators  suggests  the  possibility  of  using  such  a 
mechanism  for  subharmonic  generation.  If  we  could  find  a 
way  to  supply  power  to  the  cu-  signal  to  keep  amplitude  A, 

J  S 

constant,  the  condition  to  excite  and  signal  would  be 
governed  by  x 


dAj  2 

die’  +  >Ul>lAl  "  ff0JlA2A3  ~  °’ 
dA  <3-41 ) 

dr~  +/WktU2A2  “  ffCt>2A]A3  = 

where  r  is  properly  adjusted  to  be™ /2, 

Then  the  threshold  of  a  stable  subharrnonic  pair  with 
angular  frequencies  and  cOg  is  determined  by  t 

A~  =>Vfc^  (3-42) 

J  d 
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or  in  original  ncndimensionalized  notations  i 


2/xj/coja}^ 

L  '-n  - 


A3  "  effin 


The  correspor*  .‘ng  pressure  will  be  » 


(3-^3) 


2ji*/a)pi 


<3-V>) 


For  water,  - ^  =  6,5  J  and ,16=  0.01  poise,  the  minimum 

c0 

pressure  threshold  to  excite  the  subharnonic  pair  with  a 
1  MHZ  signal,  according  to  Eq.  (3-^*0 r  should  be  0,0017  bar. 
However,  there  are  two  factors  to  prevent  subharmonics  from 
occurring  in  real  physical  situations*  the  formation  of  shock 
wave  and  energy  spreading  into  free  space.  The  former  is 
also  caused  by  convection  and  the  nonlinearity  of  the  medium, 
as  well  as  the  accumulated  effect  of  superharmonic  generation. 
The  latter  is  due  to  the  fact  that  since  the  size  of  the 
energy  source  is  finite  and  in  real  situations  is,  of  course, 
three  dimensional,  the  wave  will  spread  out  in  part  as  a 
spherical  wave.  Under  those  conditions,  energy  is  dispersed 
faster  than  it  is  converted  to  subharmonis,  Thus,  the  phe¬ 
nomenon  of  subharmonic  generation  will  not  be  observed.  To 
avoid  such  limitations,  we  will  analyze  in  the  next  section 
the  situation  for  an  intense  wave  in  a  closed  boundary  which 


I 
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confines  the  energy  spreading.  Because  only  certain  inodes 
will  be  allowable  in  a  closed  bounaarv,  we  can  al30  eliminate 
superhamonif  generation. 


4.  Mode-Oouoling  in  a  Simple  Resonator 

In  a  closed  homogeneous  region  with  energy  reflection, 
there  exist  resonance  modes  which  arc  easily  excited.  We  will 
now  analyze  the  mode-coupling  induced  by  the  nonlinearity  of 
the  medium  in  such  a  closed  region. 

If  the  dimension  of  the  bounded  region  is  small  (say, 
about  100  wavelengths)  and  the  dissipation  and  the  nonlinearity 
of  the  medium  do  not  make  large  changes  in  wave  energy  in  a 
single  pass.  Then  we  can  consider  any  disturbance,  however 
distributed,  as  causing  the  build-up  of  a  standing  wave,  in 
this  closed  bounded  region.  The  amplitude  and  the  phase  of 
such  standing  waves  will  change  gradually  due  to  the  effect 
of  the  nonlinearity  of  the  medium*.  To  examine  their  be¬ 
havior,  we  rewrite  the  two-variable  wave  equation  derived  in 
Section  1  by  considering  slow  time  as  the  only  slow  variable. 


32u, 


ax; 


a2un 

X 

9 


=  0, 


a2u. 


3x 


32u 


3t 


9t 


2  ^1U3 


rPiui+ 


ax?  11  3tf 


ux) 
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r  B2  e2  9 2  32 

+  +  U1  **  ixf3t/l 


-A — 2 - ut  =  0  . 

3Xpt+,  * 


General  solutions  of  Eq.  (3-45)  are  s 


(3-46: 


U1  =  ^  V  V  co=  [“Vf  *  Vf  *  4(ts>j 
+  W  cos  [“Vf  +  Vf  +  W]  |  • 


where  by  our  previous  nondimensionalization, 

ol  =  k 
n  n 

If  u^  vanishes  at  boundary  xf  =  o,  then 

AA<V  =  "  Bn(ts>  ' 

W  =  W  • 


(3-4?) 


(3-48) 


(3-49) 


Eq.  (3-47)  will  be  in  a  form  of  a  standing  wave  t 

oo  r  «< 

U1  =  B=1  W  sin  [“Vf  +  ^•‘s1  J  sin  Vf  * 

(3-50) 

The  particular  values  of  kn  are  determined  by  the  end  condition 
at  x~=L?  thus,  if  u  vanishes  at  x„=L. then  k  =n^/i*»  If  the 

j  in 

end  termination  is  taken  to  be  some  frequency  dependent 
complex  impedance,  the  modal  frequencies  will  not  be  in- 


*The  effect  of  moving  boundary  is  discussed  in  Appendix  A. 
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tegrally  related* 

■'*'e  now  suppose  that.  in  a  particular  case,  there  are 
only  three  resonant  modes  of  angular  frequencies  to 2»  and 
LOj,  which  have  significant  amplitude  and  that  these  modes 
obey  the  forgoing  sum  rule  * 

a*!  +oj2  =ojj  .  (3-51) 

Considering  these  basic  modes,  we  can  rewrite  solutions  for 
Uj  in  the  form  t 

U1  -  =  ,  V'V  sin  [' "»‘f  +  'fr«<’fcs)]  3iR  kr,xf  > 

n~±  *- 

(3 ~52) 

with  the  corresponding  density  variation* 

*1  “  *  .  W  003  (“Vf  +  003  Vf  , 

n=i  «-  J 

(3—53 ) 

We  use  the  same  argument  discussed  in  Section  2  to  find  the 
second  order  function  from  Eq„  (3-^6)*  By  setting  the 
secular  terms  equal  to  zero,  we  obtain  the  results  * 


9ts 

+  jW^A-1  ■ 

‘  ffaiA2A3 

sin  r  - 

0  , 

O 

ats 

+  ^2A2  * 

•  cfc^A  1^3 

sin  r  - 

0  , 

dA- 

2 

+  jjUt^Ay  • 

sin  r  = 

0  , 

1 


t  r 
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Aiir  -  C*iA2A3  cos  r  =  0  • 
3^2 

A29t^  “  £J&;2A1A3  cos  r  =  0  , 


(3-54) 


3/^0 

A35t;  ~  ^AlA2  cos  r  =  °» 


where 


0  =  |  (  1  +  D  , 


r  =  /f,  -  /,  -  / 


2  * 


Since  the  A*s  and  the  /* s  are  functions  of  slow  time  t  only, 

s 

those  partial  differential  operators  can  be  replaced  by  total 
differential  operators.  After  combining  the  last  three  equa¬ 
tions  in  Eq,  (3-54),  we  have  mode -coupling  equations  for  a 
closed  region  (resonator). 


dA, 

2. 

«; 

+ >"iAi  - 

ao^AgA-j  sir 

r  -  0, 

dA2 

dts 

- 

ffo^A^A^  sin 

r  =  0, 

dts 

+>i^A3  + 

crcc^A^Ag  sin 

r  =  0, 

dr 

A, 

+  (  — JL  .4  .  j 

2  c^1A2A3 

^2AlA3 

dts 

A3 

A1  “ 

A2 

(3-55) 


Eq*  (3-55)  shows  that  mode -coupling  in  a  resonator  is  similar 
to  that  of  oscillators  described  in  Chapter  II  provided  the 
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coupling  parameters  in  Eq.  (2-12)  are  made  equal  in  all 
equations.  Since  the  oscillators  are  resonance  modes  of  the 
same  medium,  as  in  the  resonator,  the  result  is  less  general 
than  that  of  separate  mechanical  oscillators  which  may  have 
different  impedances  and  different  dissipation. 

For  a  nonviscous  medium,  we  can  obtain,  analogous  to 
Chapter  II,  the  three  constraints  for  the  resonator  r 

2  ?  2 
A1  +  A2  +  A3  =  e* 


cos  r  =  K. 

Since  the  system  is  assumed  to  be  lossless,  once  the 
disturbance  i?  excited,  the  energy  will  keep  exchanging 
among  the  modes.  If  we  take  the  dissipation  into  consider- 

p 

at ion,  the  teems  tt>nAn  in  Eq,  (3-55)  will  cause  the  distur¬ 
bances  to  decay  and  so,  also,  the  exchange  rate  among  the 
modes.  Eventually,  the  response  will  vanish. 

The  constraints  of  Eq.  (3-5&)  -an  be  used  to  illustrate 
the  resoonse  of  a  lossless  resonator,  V/e  do  this  by  using  a 
three  dimension  phase  space  shown  in  Fig.  3-1  which  is  similar 
to  Fig.  2-2  but  with  A1#  Ag,  and  A^  as  axes.  Since  all  the 
modes  have  the  same  impedance,  the  motion  is  eonfined  to  a 
spherical  surface  while  in  Fig.  2-2  the  surface  is  ellipsoidal. 
The  locus  in  the  pnase  space  suggests  that,  in  a  resonator, 


FIG.  3-1  PHASE  SPACE  FOR  MODE-COUPLING  IN  A 
RESONATOR 
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the  high  frequency  mode  is  quite  unstable  and  its  energy 
can  be  easily  converted  to  low  frequency  inodes.  With  dis¬ 
sipation  elements  in  the  resonator,  there  exists  a  thresh¬ 
old  for  exciting  such  subharmonic  components. 

In  general,  there  are,  in  addition  to  viscous  losses, 
loss  mechanisms  such  as  dissipation  at  the  reflecting  ends  or 
radiation  losses  along  the  sides  of  the  active  interferometer 
volume.  As  long  as  these  losses  are  sufficiently  small,  they 

may  be  considered  equivalent  to  an  effective  loss.  Thus  in 

o 

Eq.  (3-55).  we  replace  >t^nAn  by  »^nAn.  These  coefficients 
can  be  determined  directly  by  small  amplitude  measurements 
of  modal  decay  rate,  since  under  such  condition 


dA 

dt~  7nAn 


(3-5?) 


Then,  for  subharmonics  to  be  excited,  we  requires 


and 


iffi 

Aldts 


X  dA2 
A2dts 


»  0, 


or 


h 


ApAo 

<  n-fr 


A,  A- 


>»0,  or  ?2  <  , 


(3-53) 

(3-59) 


with  3in  r  =  1. 

We.  therefore,  arrive  at  a  condition  of  the  minimum  thresh¬ 
old  for  subharmonic  excitation* 


_  niHi 

3  3  * 


£q,  (3-60)  indicates  that  the  threshold  is  inversely  proper- 
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volume.  As  long  as  these  losses  are  sufficiently  small,  they 
may  be  considered  equivalent  to  an  effective  loss.  Thus  in 
Ea.  (3-55),  we  replace  >V»JnAn  by  i^An.  These  coefficients 
can  be  determined  directly  by  small  amplitude  measurements 
of  modal  decay  rate,  since  under  such  condition 


dA 


n 


dtg  "  7nAn 


A.  .  (3-57) 

Then,  for  subharmonics  to  be  excited,  we  require* 

(3-53) 


i  dAi 


Ajdt;  °-  °r  h  < 


and  1  dA2 


A,  A- 


°-  or  ?2  ■ 


(3-59) 


with  sin  r  =  1. 


We,  therefore,  arrive  at  a  condition  of  the  minimum  thresh¬ 


old  for  subharmonic  excitation* 


a3  = 


5 


(3-60) 


Eq,  (3-60)  indicates  that  the  threshold  is  inversely  propor- 
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tional  to  the  square  root  of  the  product  of  the  subharmonic 
paired  frequencies.  Hence,  when  the  paired  subharmonic 
frequencies  are  equal,  the  required  threshold  will  be  mini¬ 
mum.  Therefore  subharmonics  of  one-half  can  be  observed 
relatively  easier  provided  that  the  one-half  frequency  is 
indeed  close  to  the  resonance  mode  of  the  system. 

5.  Modes  of  a  Composite  Resonator  and  the  Selection  Rule 
for  Subharmonic  Generation 

The  results  of  the  previous  analysis  indicate  that  a 
resonator  containing  a  nonlinear  medium  can  provide  a  mecha¬ 
nism  for  subharmonic  generation.  We  have  shown  that  if  a 
resonator  has  three  resonance  modescoj-,  u>|,  coj,  and  ini¬ 
tially  there  is  a  strong  signal  in  the  mode  with  angular 
frequency  a )y  then  seme  transient  disturbance  may  cause 
other  modes  oo£  or  to£  to  appear  suddenly.  Such  a  disturbance 
may  be  weak,  but  the  energy  in  mode  ooj  will  provide  a 
source  for  modes  co*  and  u)’  through  the  nonlinear  coupling 
effect  provided  the  condition  is  approximately 

met. 

For  a  simple  one-dimensional  resonator  of  Length  L 
with  the  boundary  conditions  p1  =  0  or  u*  =  0  at  both  ends, 
the  angular  frequency  of  the  fundamental  mode  is 


(3-61) 
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ami  any  integer  multiple  of  u)^  will  still  be  a  resonance 
mode  of  such  a  system.  Therefore,  there  are  many  possibili¬ 
ties  for  modes  that  satisfy  the  condition  +  to£. 

However,  in  practical  situations,  some  means  must  be  pro¬ 
vided  for  a  driving  source.  The  actual  construction  of  a 
resonator  is  very  complicated.  Pig.  3-2  illustrates  one 
simplified  arrangement. 

To  find  the  resonance  modes  of  this  composite  reso¬ 
nator,  we  not  only  have  to  know  the  boundary  conditions  on 
both  ends  but  we  also  have  to  match  the  acoustic  impedance 
at  the  interface  due  to  two  different  media  to  provide  the 
continuity  of  pressure  and  particle  velocity.  Such  re¬ 
quirements  result  in  a  characteristic  equation  for  modes  of 

* 

this  composite  resonator  t 

R(tan  k'D  +  tan  k*H)  +  tan  k‘S 

V e  W 

-  R2  tan  k'D  tan  k'H  tan  k'S  =  0  (3-62) 

w  w 

where  R  is  the  ratio  of  wall  impedance  to  medium  impedance, 
the  k's  are  wave  numbers  in  the  corresponding  region,  and 
D,  H,  and  S  are  the  thickness  of  the  two  end  walls  and  the 
length  of  the  actual  resonator,  respectively.  Eq.  (3-62) 
can  be  rewritten  in  nondimens ional  form: 


*  Eq,  (3-62)  is  derived  in  detail  in  Appendix  B 
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FIG.  3-2  A  THREE- SECTION  COMPOSITE  RESONATOR 
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+  tan  y%  v 
s 

*r)(tan  £7  -rrHtan  |nr  ti)  =  0,  (3-68) 

h  s 


correspond  to  the  fundamental  resonance  frequency  of  each 
individual  region,  respectively.  Exact  analytical  solution 
for  the  roots  in  Eq,  (3-63)  are  impossible  to  obtain.  How¬ 
ever,  such  a  transcendantal  function  does  show  that  the 
modes  of  this  type  resonator  are  not  necessarily  integrally 
related,  as  in  the  case  of  the  simple  resonator.  Hence, 
there  may  be  only  a  limited  number  of  modes  that  will  meet 
the  (0^  =  coj|  +  u3|  condition. 

After  the  modes  of  a  resonator  are  determined  either 
by  numerical  methods  or  from  experimental  measurements,  a 
graphical  method  can  be  used  as  selection  rule  for  modes 
which  can  be  excited  as  subharmonics.  Fig,  3-3  shows  oo’ 

andCO^  used  as  vertical  and  horizontal  axes  respectively. 
The  vertical  and  the  horizontal  parallel  lines  are  drawn 
from  resonance  frequency  points  on  the  and  axes.  An 
inclined  line  connects  the  point  on  the  vertical  axis  at 
and  the  point  on  the  horizontal  axis  at  =t 


R(tan  77  or 
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GENERATION  MODES  BY  GRAPHICAL 
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Any  intersection  point  of  three  lines  gives  the  modes  with 
the  condition  =••  pj*  +ooi,.  Fig  3-3a  is  the  case  for  a 

simple  resonator.  It  indicates  that  all  the  modes  have  the 
integer  relation  and  there  are  many  intersection  points  which 
will  satify  the  condition  ^i§*  3-3b  -is  an 

example  for  the  composite  resonator  described  in  Fig.  3-2. 

As  the  figure  shows,  there  are  only  limited  number  of  modes 
meeting  the  condition  +  u>|. 


This  graphical  method  can  be  extended  to  more  compli¬ 
cated  situations.  Modes  of  a  given  resonator  can  be  de¬ 
termined  either  through  numerical  solution  or  by  actual 
experimental  measurement.  This  information  can  be  plott"'’’ 
rm  a  graph.  The  closeness  of  intersection  will  be  one  of 
the  important  factors  determining  the  possibility  of  gener¬ 
ating  subharmonics  in  the  system.  Any  intersection  point 
close  to  the  inclined  line  may  also  be  a  mode  that  can  be 
excited  if  the  driving  source  is  intense  enough  (the  detuning 
problem  will  be  further  discussed  in  next  section).  Another 
effect,  associated  with  dissipation,  is  that  there  will  be 
some  dispersion.  That  is,  at  high  frequencies  the  sound 
speed  in  the  medium  will  vary  with  frequency.  Condition 


CO.*  =  oo|  +  uJ£  docs  not  necessarily  mean  k^  -  +  k2> 


This  makes  it  more  difficult  to  excite  high  frequency  sub 


harmonics 
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6.  Detuning  for  Forcing  Excitation  of  Subharmonics  in  a 
Resonator 

If  three  modes  with  angular  frequencies  co^,  ojg,  and 
(O3  in  a  resonator  do  not  exactly  match  according  to  the  sum 
rule  of  Eq.  (3-51)»  we  then  define  a  detunig  factors  's  as? 

^01  =  ^  4  *6l' 

OJq2  ~  ^  ^  6 2» 

c0q3  =  ^3  +  €•  63,  (3-64) 

^3  /  «>x  ♦  W2, 

but  OJ03  =  «O01  4  W02. 

The  solution  of  the  first  order  wave  equation  of  Eq.  (3-45) 
is  assumed  to  be: 

u  =  nil  W  sin[CJ0ntf  +  W]  sin  k0nxf  (3'65) 
for  the  standing  wave  in  the  resonator. 

The  A's  and  fi's  are  functions  of  slow  variable  t  and  can  be 

s 

determined  by  eliminating  the  secular  terms  in  the  second 
order  equation  of  Eq,  (3-46),  The  following  is  the  result 
for  a  dissipationless  resonator: 

dA- 

dt  4  OAlA2%3  sin  r  =  °» 
s 


at :  -  ak2h;%i  sin  r  =0- 
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where  all  notations  are  defined  as  previously, 

Eq ,  (3-66)  indicates  that  the  detuning  factors  6’s  only 
appear  in  the  phase  equation.  According  to  the  analysis  in 
Section  2,  Chapter  II,  the  effect  of  the  detuning  factors 
6's  is  to  impose  a  limitation  on  the  range  of  energy  ex¬ 
change  among  the  modes. 

For  a  dissipative  resonator,  we  have  to  provide  an 
external  source  for  compensating  the  losses.  Suppose  such 
a  source  is  supplying  energy  at  one  end  of  the  resonator 
by  displacing  the  boundary  at  x  =  0  with  a  particle  velocity 
of  cos  (oj^t).  As  the  standing  wave  in  a  resonator  is 


the  combination  of  forward  and  backward  travelling  waves, 
the  propagation  of  the  particle  velocity  Uq  by  reflecting 
back  and  forth  inside  the  resonator  will  build  up  the  amp! i~ 
tude  of  a  standing  wave.  If  the  frequency  of  the  driving 
source  happens  to  be  one  of  the  resonance  modes  of  the 
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resonator,  for  each  round  trip  the  travelling  wave  propa¬ 
gates,  the  amplitude  of  the  standing  wave  will  increase  by 
2Uq.  In  mathematical  expression,  we  haves 


AA‘  2Uq 
AT*’’  =  ,  L* 


(3-6?) 


where  L*  is  the  length  of  the  resonator,  c^  is  the  sound 
speed  and.  L'/Cg  corresponds  to  the  time  taken  for  a  wave 
travelling  across  the  resonator. 

If  the  frequency  of  the  driving  source  is  only  near 
the  resonance  mode  of  the  resonator,  then  there  will  be  some 
phase  shift  for  the  standing  wave  during  the  building  up 
process.  We  can  consider  that  the  travelling  wave  contains 
in-phase  and  out-phase  components  with  respect  to  the  origi¬ 
nal  standing  wave.  Fig.  3-4  is  a  vector  diagram  showing 
such  relations.  The  in-phase  component  increases  the  ampli¬ 
tude  for  the  resulting  standing  wave  while  out-phase  compo¬ 
nent  shifts  their  phase  relation.  That  is,  for  each  round 
trip. 


A  A’  2Uq 
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Fig.  3-^  Build  up  the  Standing 
Wave  in  the  Resonator 

By  changing  to  the  nondimens ional  notation  in  slow 
rime  scale,  we  can  incorporate  the  relation  in  Eq.  (3-63) 
into  Eq.  (3-66)  to  describe  the  response  of  the  forcing 
excitation  of  subharmonic  modes  in  a  dissipative  resorator. 
The  results  are* 

dA3 

+  7qA3  +  oaia^°o3  s*n  r  ”  ^0  cos  ^3  = 


dAt 

+  YjAj  -  sin  r  =  0, 

dA2 

dt“  +  ^2A2  ~  OA1A3W02  sin  r  =  °» 

c  ■' 


A^  +  6^3  +  c?A^ AgCA^Q^  cos  r  +  GUQ  sin  /Sj  =  0, 


A.  ^  +  «1A1  ’■  ’V^Ol  005  r  =  °> 


I 


3-35 


Vfhere 


d^2 

Ag  +  S^A^  cfA^A^(i>A9  cos  r  =  0, 


TTT02 


wrc6 


(3-69) 


For  steady  state,  when  the  intensity  of  waves  in  the  resonator 
is  not  strong  enough  to  cause  the  nonlinear  effect  (a  =  0), 
the  amplitude  of  the  standing  wave  in  the  resonator  can  be 
deduced  from  Eq.  (3-69), 


A3  = 


A1  =  0, 

Ag  ~  °* 


(3-?0) 


For  a  strong  external  exitation,  the  coupling  caused  by  the 
nonlinearity  of  medium  will  result  in  a  different  set  of 
steady-state  amplitudes  of  modes  in  the  resonator* 

.2  _  (7l72  +  62  6  2) 
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A  =  —  = 

H  = 

61  *62 

h 

7  2 

7i  *7? 

Since  amplitudes  are  real  quantities,  a  minimum  value  of  UQ 
is  required  for  A^  and  Ag  to  appear.  That  threshold  for 
subharmonic  generation  is  ^iven  byt 


Wl\  *  If z2 

a2°2<O01%2  " 


(3-72) 


As  the  threshold  is  inversely  proportional  to  the  square 
root  of  the  product  of  subharmonic  frequency  pairs,  it  will 
be  a  minimum  when  they  are  equal  (U)^  =  co02). 

In  general,  only  the  acoustic  pressure  is  accessible 
for  actual  measurement.  If  component  is  taken  as  a 
reference  fcr  monitoring  the  occurrence  of  subharmonics,  we 
can  modify  A^  in  Eq.  (3-71)  in  term  of  pressure  and  express 
it  as  a  function  of  the  original  parameters  as: 


p 

r3TH 


— ^Ql,0,_  {Jk  _I  *  11.  &_2„)£ 

PVO  \  *2  4)1  4)2 

1  -.2 


(3-73) 


where  Q^,  and  Q2  are  the  quality  factor  associated  with  modes 
U)'r  and  respectively. 
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Chapter  IV 

EXPERIMENTAL  INVESTIGATION  ON  SUBHARMONIC  GENERATION 
IN  AN  ACOUSTIC  RESONATOR 

In  this  chapter,  we  describe  the  experimental  set  up, 
procedure  and  results  of  an  investigation  on  subharmonic 
generation  in  an  acoustic  resonator. 

The  acoustic  resonator  used  in  this  experiment  is  con¬ 
structed  in  the  form  of  an  interferometer  with  two  quartz 
crystals  as  the  resonator  boundaries.  One  of  the  crystals 
is  also  used  as  a  transmitting  transducer  to  supply  the 
driving  signal  into  the  resonator  and  the  other  serves  as  a 
detector  in  monitoring  the  response  of  the  system,  A  call** 
brated  acoustic  probe  smaller  Than  the  acoustic  wavelength 
of  the  signal  in  the  medium  has  been  employed  for  absolute 
acoustic  pressure  measurements  during  the  experimental  inves¬ 
tigation. 

Throughout  the  entire  experiment,  distilled  water  has 
servsd  as  a  medium.  It  has  been  filtered  and  degassed 
before  use.  The  whole  acoustic  resonator  is  enclosed  in  a 
sealed  container  so  that  contamination  of  water  during  the 
experimental  process  can  be  reduced  to  a  minimum. 

The  experimental  system  can  be  operated  in  The  frequen¬ 
cy  range  from  1  MHZ  to  5  MHZ,  although  detailed  studies  in 
observing  subharmonics  are  only  carried  out  around  the  re- 


gion  oi'  1.5  MHZ,  which  is  the  natural  resonance  frequency  of 
the  transmitting  quartz  crystal.  As  the  validity  of  using 
the  acoustic  probe  is  limited  in  measuring  of  absolute 
acoustic  pressures  at  high  frequency,  the  subharmonic  pheno¬ 
menon  is  only  considered  for  qualitative  reference  when  the 
crystal  is  driven  around  4,5  MHZ,  the  third  harmonic  of  its 
natural  frequency. 

Data  concerned  with  actual  modes  of  the  resonator,  the 
loss  factor  associated  with  each  mode,  and  the  signal  thresh¬ 
old  for  exciting  subharmonics  at  various  resonator  lengths 
and  at  different  signal  frequencies  are  recorded.  Such  ^for¬ 
mation  has  been  compared  with  the  theoretical  analysis  dis¬ 
cussed  in  the  previous  chapters. 

1,  Conditions  for  the  Experimental  Investigation  of  Sub- 
harmonie  Generation 

According  to  the  previous  analysis,  there  are  two 
necessary  conditions  for  subharmonics  to  be  excited.  One 
is  that  the  system  should  contain  the  subharmonic  modes  and 
the  other  is  that  the  high  intensity  acoustic  wave  should  be 
built  up  to  provide  the  necessary  nonlinear  coupling  mech¬ 
anism. 

With  this  consideration,  an  interferometer  type  reso¬ 


nator  is  adopted  for  our  experimental  investigation  of  sub¬ 
harmonic  generation.  In  general,  the  interferometer  pos- 
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sesses  many  resonance  inodes.  If  one  such  resonance  mode  is 
excited  by  a  external  source,  a  high  intensity  standing  wave 
can  easily  be  achieved.  Any  modes  whose  frequencies  are 
below  the  external  driving  signal  can  ther  be  considered  sub¬ 
harmonic  modes.  Another  advantage  for  using  an  interfer¬ 
ometer  as  a  resonator  is  that  the  end  reflectors  can  be  util¬ 
ized  as  transducers  to  provide  a  means  of  supplying  the  driv¬ 
ing  signal  or  detecting  the  system's  response.  Since  there 
is  no  side  boundary  to  enclose  the  interferometer,  only  one 
type  of  mode  exists  in  such  a  resonator.  In  the  case  that 
the  wavelength  of  the  acoustic  v/ave  in  the  medium  is  much 
smaller  than  the  dimension  of  the  end  reflectors,  the  inter¬ 
ferometer  v.  ill  act  almost  as  a  one-dimensional  resonator. 

The  design  of  interferometer  depends  very  much  on  its 
operating  frequency.  Most  interferometers  used  for  experi¬ 
mental  observation  of  subharmonics  are  in  the  frequency 
range  3  MHZ  to  5  -9,38^  g-he  advantage  of  apply¬ 

ing  high  frequency  signal  for  subharmonic  generation  is  that 
the  threshold  for  exciting  subharmonics  is  lov/er  according 
to  the  inverse  relation  between  the  signal  frequency  and 
the  threshold  derived  in  our  previous  analysis,  Eq,  (3~6Q) 
Chapter  III .  An  additional  reason  for  working  in  this  fre¬ 
quency  range  is  that  a  ham  radio  transmitter  can  be  easily 
modified  for  the  driving  signal  power  source.  But  there  are 
some  disadvantages  to  using  high  frequency  signals.  The 
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important  one  is  that  the  absolute  instantaneous  acoustic 
pressure  will  be  difficult  to  measure  without  causing  scat¬ 
tering  effects.  (The  wavelength  in  that  frequency  range  will 
b°  less  than  0.03  cm  in  water.)  The  laser  beam  diffraction 
method,  employed  by  previous  experimenters  as  a  ;ans  of 
detecting  the  acoustic  wave,  does  not  directly  give  the 
absolute  acoustic  pressure  and  it  involves  elaborate  equip¬ 
ment  and  alignment.  A  second  handicap  for  driving  the  inter¬ 
ferometer  at  high  frequencies  is  that  a  thinner  quartz 
crystal  plate  has  to  be  used  for  the  transmitting  transducer. 
This  presents  the  difficulty  of  mounting  the  crystal  for 
direct,  coupling  with  the  medium.  A  common  practice  is  to 
cement  the  crystal  on  a  metal  plate  for  making  a  composite 
transducer-^.  Our  preliminary  tests  of  such  an  arrangement 
shows  that  when  an  intense  signal  i.s  applied  to  the  compos¬ 
ite  transducer,  the  cement  between  crystal  and  metal  plate 


seems  to  change  its  adhesive  force,  and  a  nonlinear  response 
begins  to  appear.  In  some  cases,  even  subharmonics  are  pre¬ 
sented.  They  are  very  stable  and  can  not  be  affected  by  .just, 
changing  the  length  of  the  interferometer. 


Surveying  the  current  literatures,  we  found  that  it  is 

4  o, 

possible  to  make  an  acoustic  probe  of  size  about  0.0?.  cm 

ki 

,  This  corresponds  to  the  wavelength  of  a  5  MHZ  acoustic 


wave  in  water.  Then  the  proper  highest  frequency  that  should 
be  used  without  causing  seriously  scattering  effects  will  be 
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around  1.5  MHZ,  which  requires  a  0.19  cm  thick  quartz  plate 
for  the  transmitting  transducer  in  the  interferometer. 

Quartz  of  such  thickness  is  strong  enough  to  be  used  as  the 
reflecting  wall  in  the  interferometer.  This  solves  the  prob¬ 
lem  by  avoiding  the  complication  of  using  a  comocsite  trans¬ 
ducer,  A  much  lower  frequency  quartz  is  chosen  for  the 
other  reflecting  wall  of  the  interferometer  for  two  reasons i 
first,  it  will  have  better  sensitivity  in  detecting  the  sub¬ 
harmonic  components;  second,  its  thickness  will  allow  for  a 
knife  edge  mounting  for  the  crystal,  thereby  reducing  me¬ 
chanic  losses. 

In  our  first  pilot  interferometer  composed  of  1.5  MHZ 

and  600  KHZ  quartz  crystal  plates  as  the  reflecting  walls  of 

the  interferometer,  the  Q  (quality  factor)  of  such  a  reso- 

3 

nator  measured  during  preliminary  test  was  about  10  .  Ac¬ 
cording  to  the  relation  derived  in  Chapter  III,  Eq.  (3-73), 
an  acoustic  pressure  amplitude  of  at  least  15  bars  is  needed 
for  exciting  subharmonics  for  this  Q. 

An  acoustic  standing  wave  with  this  intensity  will  also 

42 

cause  cavitation  in  ordinary  water  .  To  avoid  the  compli¬ 
cation  of  the  cavitation  phenomenon  during  the  observation 
of  subharmonics,  we  decided  to  process  the  liquid  before 
using  it  in  the  interferometer.  The  threshold  for  cavitation 
will  increase  after  the  liquid  has  been  purified  through  a 
degassing  and  filtering  process.  Hence,  the  subharmonic 
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threshold  will  be  less  than  the  cavitation  threshold. 

We  also  designed  a  system  to  enclose  the  interferometer 
so  that  the  treated  liquid  would  not  be  subsequently  contami¬ 
nated,  This  involves  the  elaborate  construction  of  a  sealed 
system.  As  we  do  not  find  that  previous  experimenters  speci¬ 
fied  how  they  have  treated  the  liquid  used  in  their  investi¬ 
gation  of  subharmonics,  we  hoped  that  our  effort  to  make  the 
whole  test  system  a  closed  one  would  provide  some  new  and 
interesting  information. 

There  are  other  advantages  to  using  a  low  driving 
signal  frequency  in  the  resonator.  Since  the  wavelength  of 
acoustic  wave  in  the  water  is  longer,  the  mechanical  align¬ 
ment  for  parallelism  of  the  two  reflecting  plates  is  easier 
and  there  are  less  appreciable  detuning  drifts  due  to  temper¬ 
ature  changes  in  the  interferometer. 

Based  on  the  above  considerations,  our  final  design 
of  the  test  system  for  observing  subharmonic  generation 
phenomena  possesses  the  following  special  features* 

(a)  The  acoustic  transducers  which  also  serve  as 
reflecting  walls  for  the  interferometer  are  directly  coupled 
to  the  medium  in  the  resonator, 

(b)  The  medium  used  in  the  resonator  has  been  care¬ 
fully  filtered  and  degassed  to  a  high  purity.  Cavitation  is 
avoided  in  the  range  of  tests  used  for  observing  the  subbar- 
monic  phenomenon. 
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(c)  An  acoustic  probe  smaller  than  half  a  wavelength 
of  the  acoustic  wave  in  water  can  be  calibrated  for  record¬ 
ing  the  instanxaneous  acoustic  pressure. 

(d)  The  system  can  be  operated  in  a  pulsed  mode 
for  studying  the  transient  behavior  of  subharmonic  genera¬ 
tion. 

Pig.  4-1  is  the  block  diagram  of  whole  testing  system. 
In  next  section,  each  component  will  be  described  in  detail. 


Matching 

Network 


t 


FIG.  4-1  BLOCK 


!\GRAM  FOR  EXPERIMENTAL  INVESTIGATION 
OF  SUBHARMONICS 


4-9 


2 .  Apparatus 

The  actual  apparatus  used  in  this  experiment  for  the 
investigation  of  subharmonic  generation  can  be  seen  in  the 
photograph  of  Pig.  4-2.  According  to  function,  we  can 
divide  the  apparatus  into  three  subsystems!  resonator,  the 
liquid  purification  equipment,  and  electronic  instruments. 
Some  items  are  a^'jpted  from  available  commercial  products 
but  most  equipment  is  designed  and  constructed  only  for  our 
special  purpose  in  this  experiment.  Each  subsystem's  func¬ 
tion  and  construction  is  described  belowt 

(a)  The  Resonator 


The  construction  of  the  resonator  is  shown  by  the 
sketch  in  Fig.  4-3.  It  mainly  contains  an  interferometer 
with  two  x-cut  quartz  crystal  plates  as  its  boundary  walls. 
The  space  between  these  two  walls  can  be  adjusted  from  0  to 
13  cm.  The  fine  adjustment  is  controlled  by  a  micrometer 
which  only  can  cover  a  range  of  2,5000  cm.  An  additional 
coarse  adjustment  is  provided  for  larger  space  changes. 

One  of  the  quartz  plates,  which  we  use  as  a  trans¬ 
mitting  transducer  to  supply  the  acoustic  signal" to  the 
resonator,  is  a  disc  of  2.54  cm  in  diameter  and  0.19  cm  in 


thickness  whose  natural  resonance  frequency  in  air  is  1.5 


MHZ.  One  side  of  this  disc  is  plated  with  gold  on  chrome  and 


the  other  side  is  plated  with  the  same  material  but  is  left 
with  a  0,0625  cm  wide  blank  annulus  on  the  edge.  This  disc 
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FIG.  4-3  THE  RESONATOR 
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is  mounted  on  a  brass  housing  with  an  open  hole  of  2,2  cm  in 
diameter  so  that  the  large  part  of  the  plated  surface  will 
be  in  direct  contact  with  the  liquid  in  the  resonator.  A 
plexglass  ring  is  used  on  the  other  side  of  the  quartz  disc 
for  pressing  it  against  the  brass  housing.  Silicon  grease 
is  applied  around  the  edge  of  the  quartz  disc  before  instal¬ 
lation  so  that  there  is  a  good  seal  between  disc  and  brass 
housing  to  prevent  a  leak.  Electrical  contact  with  the  back 
face  of  the  quartz  disc  is  made  with  a  sc-ft  spring  which  is 
soldered  to  a  BNC  connector.  During  operation,  the  signal 
is  introduced  through  the  BNC  connector  and  the  brass  housing 
and  the  front  surface  of  the  quartz  disc  are  at  ground  poten¬ 
tial, 

A  much  thicker  quartz  disc  is  used  for  the  other 
boundary  wall  of  the  interferometer  and  acts  as  a  detector 
for  the  system's  response.  This  disc  also  has  a  diameter 
of  2,5^  cm  and  is  clamped  on  the  knife  edge  of  a  brass  ring 
which  can  be  mounted  on  an  another  brass  nousing.  Between 
the  knife  edge  and  the  circumference  of  the  quartz  disc, 
stycast  2850  FT  (manufactured  by  Emerson  and  Cuming,  Inc.) 
is  used  to  avoid  possible  leakage.  One  side  of  this  unit  is 
vacuum  deposited  with  silver  as  the  front  surface  of  the 
boundary  wall.  The  other  side  of  the  quartz  disc  is  plated 
with  gold  on  chrome  and  touches  a  soft  spring  which  serves 
as  one  of  the  electrodes  after  this  unit  is  mounted  on  a 
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brass  housing.  An  0  ring  is  inserted  between  the  brass 
housing  and  the  detecting  quartz  unit  for  leakage  prevention. 
We  have  made  several  such  units  with  various  natural  reso¬ 
nance  frequency  quartz  discs:  1  MHZ,  600  KHZ,  and  JoO  KHZ. 
They  are  easily  interchangeable.  The  housing  of  the  de¬ 
tecting  unit  is  a  part  of  a  brass  telescope  which  has  a  mi¬ 
crometer  to  adjust  its  length. 

Both  transmitting  transducer  and  detecting  unit  are 
mounted  separately  on  2  cm  thick  brass  plates.  A  pyrex 
cylinder,  10  cm  inside  diameter  and  1 ^  cm  in  length,  with 
two  teflon  gaskets  is  inserted  between  these  two  brass 
plates  to  form  a  container.  The  whole  unit  is  fixed  on  a 
steel  frame.  There  are  three  adjustment  bolts  for  the 
alignment  of  the  parallelism  between  the  transmitting 
transducer  and  the  detecting  unit.  An  ir.let  and  an  outlet 
are  provided  for  this  container  for  filling  and  circulating 
the  liquid.  Fig.  4-4  is  the  photographical  close  view  of 
the  actual  arrangement. 

an  acoustic  probe  is  also  placed  inside  the  pyrex  con¬ 
tainer  for  absolute  acoustic  pressure  measurement.  The 
active  tip  of  this  probe  is  about  0,05  cm  in  diameter  and 
has  a  microdot  connector  on  its  end  (Microdot  Inc,,  part 
No.  032  0015  0001).  This  unit  is  then  attached  to  a  0,5  cm 
diameter  brass  tubing  going  through  the  mounting  plate  with 
an  0  ring  as  a  sealing  element.  The  probe  tip  is  detachable  1 


FIG.  4*5  THE  ACOUST.'C  PROBE 
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the  picture  in  Fig.  4- 5b  shown  its  actual  size. 

The  construction  of  the  acoustic  probe  is  shown  in  th° 
detailed  sketch  in  Pig.  4-5a«  The  tin  is  BaTiO^  rowdor 
melted  on  a  0,02  cm  diameter  pj.atinum  w ire  by  a  microflame 
torch  (manufactured  by  Microflame  Inc,).  The  main  frame  of 
the  probe  is  0.055  cm  diameter  stainless  tubing  covered  wit'"' 
a  layer  of  0.B,  Olyptal  (1201  Red  Enamel).  The  external  sur¬ 
face  of  the  probe's  main  frame  is  coated  with  B-Kote  40 
silver  paint  (manufactured  by  Epoxy  Products  Co.)  to  serve 

as  an  electrode  for  the  orobe.  The  BaTiO^  tin  is  polarized 

/ 

with  30  V  DC  at  130°C. 

The  sensitivity  of  such  a  probe  with  the  preamplifier 
is  about  23  mV  per  bar.  The  ealibration  procedure  for  the 
probe  will  be  described  later. 

(b)  The  Liquid  Purification  Equipment 
To  control  the  purity  of  the  liquid  under  experimental 
investigation,  the  container  with  its  interferometer  is  con¬ 
nected  through  Tygon  Vacuum  Tubing,  3/^"  in  inside  diameter 
(made  by  Norton  Plastics  and  Synthetics  Division ),  to  the 
accessories  for  the  liquid  cleanin'"  process.  Fig.  4-6  shows 
the  detailed  arrangement. 

Two  Mil linore  Filters  (manufactured  by  i-.ili ioore  Fil¬ 
ter  Co.)  of  5.0  m’cron  po^e  size  are  used  :n  the  cleanin'? 
system.  One  is  for  trapping  the  contaminant  before  the 
liquid  gets  into  tne  system  the  other  is  inserted  the 
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FIG-  4-6  ARRANGEMENT  AND  EQUIPMENT  FOR 
FILTERING  AND  DEGASSING 
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system  for  closed-loop  filtering. 

A  Fisher  Airejector  is  adopted  as  a  vacuum  pump  for 
degassing  the  liquid  in  the  system.  Such  a  simple  arrange¬ 
ment  can  easily  reach  a  vacuum  less  than  30  mm  Hg,  The 
vacuum  is  monitored  by  a  pressure  gauge  (made  by  Marshalltown 
MFG .  Co. ) « 

Circulation  of  the  liquid  around  the  system  is  pro¬ 
vided  by  a  micro-bellows  pump  (manufactured  by  Research 
Appliance  Co.  serial  No.  619616  4648).  During  the  actual 
experimental  investigation  of  sabharmonics,  the  operation 
of  the  pump  is  stopped  in  order  to  avoid  the  possible  inter¬ 
ference  caused  by  the  vibration. 

(c)  Electronic  Instruments 

We  have  two  signal  generators  for  use  as  the  signal 
sources  in  this  experiment.  The  General  Radio  1001-A,  which 
covers  a  range  from  5  KHZ  to  50  MHZ,  is  used  to  check  the 
resonance  modes  of  the  interferometer.  The  other,  General 
Radio  121 1-B  Unit  Oscillator,  mainly  serves  as  the  signal 
input  to  the  power  amplifier  for  the  transmitting  transducer. 
This  oscillator  can  be  modulated  by  an  external  source  through 
its  power  supply.  We  have  built  a  SOP  gate  in  conjunction 
with  a  General  Radio  Unit  Pulser  1217-A  to  modify  it  as  a 
pulsed  signal  generator.  The  transient  behavior  of  subhar¬ 
monic  generation  in  the  interferometer  can  be  studied  by 
such  conversion. 
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The  power  amplifier  was  originally  designed  and  cort- 
42 

structed  by  Barger  for  his  research  work  on  acoustic  cavi~ 
tation.  This  amplifier  can  be  driven  to  400  watts  output 
with  a  1  V  rms  input  signal  and  has  the  half  power  points  at 
12  KHZ  and  1.5  MHZ.  We  find  it  quite  suitable  for  driving 
the  quartz  transducer  of  the  interferometer  at  a  high  impe¬ 
dance  level,  A  General  Radio  Variable  Inductor  type  10?-J 
is  used  to  tune  out  the  stray  capacity  in  its  output  cir¬ 
cuit. 

During  preliminary  studies  on  subharmonic  generation 
phenomenon,  we  have  employed  a  Lysco  600-S  ham  transmitter 
as  a  power  scurce.  The  output  stage  of  this  transmitter  is 
a  beam  power  pentode  807  which  can  deliver  maximum  power  up 
to  50  watts,  A  regulated  power  supply,  Lamda  model  71,  is 
used  as  the  plate  voltage  supply  for  the  pentode  to  adjust 
the  output  level.  A  matched  network  has  to  be  inserted 
between  the  transmitter  and  the  quartz  transducer  to  obtain 
optimal  results.  This  arrangement  was  solely  for  qualitative 
studies  of  subharmonics  when  the  quartz  transducer  is  oper¬ 
ated  around  its  third  harmonic,  4,5  MHZ, 

During  the  experimental  investigation,  the  frequency 
of  the  driving  signal  is  checked  by  a  Hewlett-Parkard  52441. 
Electronic  Counter, 

For  parallelism  alignment  of  the  interferometer,  the 
signal  is  provided  by  a  Hewlett-Parkard  type  161  Pulse  Gener- 
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ator.  The  alignment  procedure  will  be  described  later. 

There  are  two  separate  channels  for  detecting  the  re¬ 
sponse  of  the  system.  The  acoustic  probe  is  mainly  for 
absolute  acoustic  pressure  measurement.  It  has  a  self  con¬ 
tained  battery  operated  preamplifier  for  eliminating  ground¬ 
ing  and  noise  problems.  Pig.  &-7  shows  its  circuitry  which 
has  about  50  dB  gain.  The  output  of  this  preamplifier  can 
be  fed  directly  to  the  scope  to  display  the  waveform. 

Another  channel  is  provided  by  the  quartz  disc  wall  in 
the  interferometer.  Since  this  detector  has  high  sensitivity 
to  the  subharmonics,  its  output  is  directly  connected  to  a 
scope,  Tektronix  type  5^5  Dual-Beam  Oscillascope,  and  a  v/ave 
analyzer,  Hewlett-Parkard  model  310A.  The  wave  analyzer 
also  serves  as  a  filter  for  picking  up  the  particular  subhar¬ 
monic  component.  The  frequency  of  this  subharmonic  is  then 
measured  by  an  another  electronic  counter,  Kewlett-Parkard 
model  523 B. 

Voltage  and  current  probes,  both  manufactured  by 
Tektronix,  are  employed  to  monitor  the  output  waveform  of 
the  power  amplifier  so  that  we  can  be  sure  that  no  distorted 
signal  drives  the  transmitting  transducer  during  experimental 
investigations. 

3.  Experimental  Procedure 

In  investigating  the  subharmonic  generation  phenomenon, 
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quantitative  measurements  can  be  made  of  two  important  pa¬ 
rameters*  the  driving  signal  threshold  for  exciting  the  sub¬ 
harmonics  and  the  frequency  components  of  the  excited  sub- 
harmonics.  These  parameters  depend  on  the  physical  property 
of  the  medium  as  we  have  concluded  in  the  previous  analysis. 
However,  in  the  actual  experimental  situation,  the  results 
are  also  effected  by  the  way  the  system  is  set  up.  We  will 
mention  some  of  our  experiences  during  preliminary  test 
in  the  early  experimental  stage  before  describing  the  actual 
procedure  we  finally  adopted  for  data  collection. 

Without  any  treatment  of  the  water  used  in  the  reso¬ 
nator,  we  found  it  was  almost  impossible  to  excite  subhar¬ 
monics  with  some  air  bubbles  attached  to  the  face  of  the 
boundary  walls  of  the  interferometer.  However,  even  after 
those  bubbles  are  removed,  either  externally  or  by  sell 
dissolution,  the  excited  subharmonics  appear  to  be  very  un¬ 
stable.  Further  increases  in  the  driving  signal  intensity 
only  seem  to  extinguish  the  existence  of  subharmonics.  This 
peculiar  phenomenon  was  finally  explained  when  we  discovered 
that  gaseous  cavitation  was  occurring  along  the  path  of  the 
acoustic  standing  wave  in  the  interferometer.  Tiny  gas 
bubbles,  whicn  can  be  observed  by  shining  collimated  light 
on  them,  are  dancing  around.  Sometimes  they  aggregate  to¬ 
gether  to  become  a  larger  bubble  and  then  move  away  by 
buoyancy  forces.  Such  activities,  we  believe,  will  absorb 
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more  acoustic  energy  and  interfere  with  the  resonance  modes 
of  the  interferometer  so  that  the  threshold  of  exciting  sub- 
harmonics  is  increased,  makes  the  sustenance  of  subharmonics 
more  difficult.  In  the  low  frequency  range,  the  threshold 
for  acoustic  cavitation  is  lower  than  the  subharmonic  thresh¬ 
old  and  therefore  the  observation  of  subharmonics  becomes 
more  difficult.  Hence  we  reached  the  conclusion  that  we 
should  control  the  water  used  in  the  experiment  tc  avoid  the 
interfering  effects  of  gaseous  cavitation. 

Because  the  energy  of  the  intense  acoustic  wave  is 
absorbed  in  the  resonator,  we  have  observed  that  subharmonics 
become  unstable  due  to  temperature  changes  during  a  long 
period  of  operation.  A  2 °C  change  in  temperature  has  be  re¬ 
corded  after  about  5  hours  continual  operation.  We  have 
overcome  this  problem  by  using  pulsed  signals  to  reduce  the 
heat  generated  in  the  water. 

Our  initial  experience  with  the  instruments  has  re¬ 
vealed  some  misleading  data  from  the  quantitative  measure¬ 
ments,  Ore  of  them  is  that  the  quartz  disc  serving  as  the 
boundary  wall  of  the  interferometer  can  only  be  used  as  a 
detector  for  the  system's  response.  This  is  a  very  conven¬ 
ient  way  to  monitor  subharmonic  phenomena  but,  as  the  re¬ 
sponse  of  the  particular  quartz  disc  is  frequency  dependent, 
it  is  not  a  proper  way  to  take  quantitative  data.  A 
similar  situation  also  aonlies  to  the  transmitting  trans- 
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ducer,  The  acoustic  intensity  inside  the  interferometer  is 
not  necessarily  proportional  to  the  voltage  across  the 
quartz  plate  but  depends  on  other  factors  such  as  alignment 
and  tuning  of  the  interferometer,  separation  of  two  boundary 
walls,  the  natural  resonance  frequency  of  the  quartz  disc, 
the  condition  of  radiative  beam  spreading,  and  losses  due  to 
the  viscosity  of  the  liquid,  etc.  This  is  true  even  when 
different  quartz  crystal  discs  are  used,  with  the  different 
driving  signal  frequencies,  since  their  thickness  and  condition 
for  mounting  are  changed.  The  mounting  loss  plays  an  im¬ 
portant  role  in  the  final  intensity  of  an  acoustic  wave  in 
the  resonator.  It  becomes  totally  improper  to  relate  a  set 
of  applied  voltages  on  the  quartz  discs  to  the  true  acoustic 
intensity  in  the  interferometer.  For  this  reason,  we  devel¬ 
oped  a  small  acoustic  probe  for  absolute  pressure  measure¬ 
ment. 

Alignment  of  the  two  boundary  walls  of  the  interfero¬ 
meter  is  very  critical  for  investigating  the  subharmonic 
phenomenon.  The  threshold  for  exciting  subharmonics  and  the 
frequency  components  of  the  excited  subharmonics  will  be 
altered  greatly  ju3t  by  manipulating  the  parallelism  of  two 
boundary  walls.  In  order  to  obtain  consistent  experimental 
results,  we  set  up  a  standard  alignment  procedure,  which 
will  be  described  later,  before  collecting  the  data. 

The  following  is  a  detailed  description  of  the  experi- 
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mental  procedure  we  have  adopted* 


(a)  Calibration  of  The  Acoustic  Probe 
The  acoustic  probe  is  primarily  designed  for  use  at 
high  frequencies.  Since  measurements  of  input  current  to 
the  transducer  at  this  frequency  range  is  not  very  accurate, 
the  reciprocity  method  for  calibration  becomes  unsuitable 
under  present  conditions.  A  simple  method  is  adopted  from 
the  current  literature  « 

In  studying  the  propagation  of  a  finite -amplitude 
wave  in  a  liquid,  one  of  the  established  results  is  that 
a  sinusoidal  wave  at  the  source  of  radiation  will  become 
a  sawxooth  at  some  distance  away  from  the  source.  This 
distance  can  be  expressed  in  terms  of  wave  parameters  and 
properties  of  the  medium  by  the  relation* 


1  = 


A  Pc 


(4-1) 


ir(Y  +  l)p  , 

where  \  is  the  wavelength,  P  is  the  density  of  the  liquid, 
c  is  the  sound  speed,  p  is  the  pressure  amplitude  of  the 
acoustic  wave,  and  y  is  the  nonlinear  parameter  of  the 
liquid  (about  6.5  for  water).  In  the  case  of  spherical 
wave  divergence  in  the  far  field  region,  the  variation  of 
the  amplitude  of  the  formed  sawtooth  wave  with  the  distance 
from  the  source  is  determined  by  the  equation: 
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v/here  f  is  the  frequency,  and  pQ  and  px  are  the  amplitudes 

of  xhe  sawtooth  at  distances  and  x  from  the  source  respec- 

o 

tively.  This  relation  can  be  rewritten  in  the  form  of  ratios 

and  x  /x  to  obtain  the  absolute  oressure  amplitude 
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Once  the  pressure  is  known,  we  can  easily  compute  the  sensi¬ 
tivity  of  the  acoustic  probe  from  its  electrical  output 
signal. 

The  calibration  of  the  probe  used  in  this  experiment 
is  carried  out  in  an  anechoic  tank  which  has  the  dimensions 


75  cm  wide,  125  cm  long,  and  75  cm  in  depth,  A  0,254  cm 
thick  PZT-4  disc:  of  2  cm  in  diameter  (manufactured  by  Clevite 
Co.)  is  used  for  the  transmitting  transducer.  The  resonance 
frequency  of  this  disc  is  about  7^0  KHZ.  For  calibrating: 
the  probe  at  high  frequencies,  we  used  the  second  harmonic 
component  of  the  distorted  wave  as  the  reference.  With  the 
preamplifier,  the  sensitivity  of  the  probe  around  1.534  KHZ 
is  23  mv/bar. 


t 
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(b)  Treatment  of  Water 

The  arrangement  and  equipment  for  the  treatment  of 
water  used  for  the  experimental  investigation  is  shown  in 
Pig.  4-6,  The  container  is  first  evacuated  partially  by 
the  Ai rejector  so  that  the  water  can  be  sucked  in  through 
the  inlet  filter.  After  the  water  has  completely  filled 
up  the  system*  we  then  shut  off  the  inlet  and  at  same  time 
close  the  system  by  proper  turning  of  a  three  way  valve. 

The  water  can  be  continually  filtered  through  the  in-line 
filter  by  the  forced  circulation  due  to  the  micro-bellows 
pump. 

The  degassing  process  is  done  after  the  water  has 
gone  through  the  filter  several  times.  A  30  mm  Hg  is 
generally  reached  during  the  final  stage  of  degassing. 

Then  we  open  up  the  relief  valve  for  the  system  and  the 
hydrostatic  pressure  in  the  container  returns  to  atmos¬ 
pheric  . 

This  process  is  repeated  at  the  beginning  of  each 
measurement  ana  at  the  beginning  of  each  day  so  we  car  be 
sure  that  no  further  contamination,  possibly  caused  by 
corrosion  and  leakage  of  the  components,  will  deteriorate 
the  system. 

During  the  experiment,  we  stop  the  circulating  pumr  to 
avoid  possible  interference  from  vibration,  A  focused  beam 
of  light  ft'orn  a  100-watt  projector  lamp  ic  used  to  aid  in 
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the  visual  detection  of  possible  gas  bubbles  in  the  inter¬ 
ferometer.  With  this  arrangement,  we  are  able  to  ob¬ 
serve  subharmonics  at  1.5  MHZ  without  the  occurrence  of  cavi¬ 
tation, 

(c)  Alignment  of  The  Interferometer 

The  arrangement  for  alignment  of  the  interferometer  is 
shov/n  in  Fig.  4-8a.  The  transmitting  transducer  is  ener¬ 
gised  by  a  10  microsecond,  50  volts  DC  pulse.  This  pulse 
will  excite  the  quartz  disc  in  the  transducer  and  an  acoustic 
wave  will  radiate  out  into  the  water  inside  the  interfero¬ 
meter.  The  pulsed  wave  will  be  picked  up  again  by  the  trans¬ 
mitting  transducer  after  it  has  reflected  back  from  the 
other  boundary  wall  of  the  interferometer.  This  reflected 
signal  can  be  detected  with  a  tuned  tank  connected  to  the 
transmitting  transducer  and  displayed  on  the  oscilloscope, 
if  the  repetition  rate  of  the  DC  pulse  is  slow  compared  to 
the  time  required  for  wave  traversal  across  the  interfero¬ 
meter,  the  successive  reflected  pulsed  signals  v/ill  also  be 
detected  by  the  scope.  By  monitoring  the  amplitudes  of  those 
reflected  pulses,  we  can  align  the  reflected  wall  of  the 
interferometer  with  its  adjusting  screws.  The  best  alignment 
for  the  parallelism  of  the  boundary  walls  is  obtained  when 
the  amplitudes  of  pulsed  signals  become  a  maximum  and  de¬ 
crease  exponentially  according  to  the  time  sequence  of  their 
arrival.  This  alignment  proceduce  is  not  affected  by  the 
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length  of  the  interferometer.  A  typical  picture  of  the 
response  signal  on  the  oscilloscope  is  shown  by  Pig,  4-Bb 
after  the  interferometer  has  been  aligned.  The  decrease 
in  the  amplitudes  of  the  successive  pulses  also  gives  the 
information  about  the  quality  factor  (Q)  of  the  resonator 
at  the  given  resonance  frequency, 

(d)  Excitation  of  Subharmonics 

A  strong  standing  wave  can  be  established  when  we 
drive  the  interferometer  at  one  of  its  resonance  modes. 

This  is  accomplished  by  adjusting  the  length  of  the  inter¬ 
ferometer  with  the  micrometer  for  a  given  input  signal  or 
changing  the  signal  frequency  for  a  fixed  length  interfero¬ 
meter.  By  increasing  the  signal  input  to  the  transmitting 
transducer  of  the  interferometer,  subharmonics  may  be  ob¬ 
served  from  the  response  of  the  detecting  quartz  disc  after 
the  intensity  of  the  acoustic  wave  inside  the  interferometer 
exceeds  a  certain  level.  A  typical  response  of  such  subhar¬ 
monic  generation  in  the  time  domain  '3  displayed  in  Pig.  4— 9a 
from  the  picture  taken  on  the  oscilloscope.  The  frequency 
components  of  such  signals  can  be  further  analyzed  through 
the  wave  analyzer?  the  exact  frequency  is  read  from  the  elec¬ 
tronic  counter.  The  meter  on  the  wave  analyzer  also  records 
the  relative  intensity  of  the  particular  subharmonic  compo¬ 
nent. 

The  frequency  range  of  the  wave  analyzer  covers  from 


f  input  =  1599=2  kHz 
f  sub  s  799.6  kHz 

(,2v/cm,  l^s/cm) 


*  f  input  s  1544.8  kHz 


fsub  s  1154.5  kHz 
390.3  kHz 
(.2v/cm,  1/iS/cm) 


(a)  SUBHARMONIC  RESPONSES 


=  390.3  kHz 

(.2  v /cm *  5ms /cm)  (.2y  /cm,  5  ms/cm) 

(b)  GROWTH  OF  SUBHARMONICS 


FIG.  4-9  OBSERVATION  OF  SUBHARMONICS 
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1  KHZ  to  1.5  MHZ.  Sometimes,  at  high  driving  signal  levels, 
more  than  one  subharmor.ic  pair  can  be  observed. 

For  gathering  information  about  the  threshold  for 
subhafmonic  generation,  we  used  long  pulse  train  of  the  driv¬ 
ing  signal  applied  to  tjhe  interferometer.  The  actual  growth 
of  the  subharmonic  acoustic  wave  can  be  displayed  on  the  os- 

1  ' 

cilloscope,  such  as  the  one  shown  on  Fig  4-9b.  From  that 
picture,  the  point  where  the  wave  starts  to  break  and  sub- 

1 

hjarmoni,cs  begin  to  appear  can  be  easily  determined.  If 
t'his  response  is  measured  using  the  output  of  the  acoustic 
pjrobe,  the  level  at  the  break  point  is  the  threshold  for 
subharmonic  excitation, 

'  (e)  Measurement  of  the  Sysrcm  Parameters  of  the 

f 

|  Interferometer 

i 

For  comparing  the  experimental  data  and  the  theoretical 

{ 

analysis,  information  concerning  the  actual  condition  of  the 
interferometer  is  examined  to  explore  its  relation  with  the 
mathematical  model  of  the  ideal  resonator.  Two  parameters 
have  been  recorded  as  the  quantitative  measurement  of  the 
physical  condition  of  the  interferometer.  The  first  parame¬ 
ter  is  the  actual  resonance  modes  of  the  interferometer  at  a 
given  experimental  condition.  The  other  is  the  dissipation 
factor  associated  with  each  resonance  mode, 

A  simple  way  has  been  adopted  for  measuring  the  reso¬ 
nance  modes  of  the  interferometer.  A  variable  signal  fre- 
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quency  generator  is  directly  connected  to  the  transmitting 
transducer  of  the  interferometer.  By  changing  the  frequency 
of  the  input  signal,  we  observe  the  amplitude  response  from 
the  detector  at  one  of  the  boundary  walls.  The  frequency 
at  which  the  response  amplitude  is  maximum  is  the  resonance 
mode  frequency. 

For  determining  the  dissipation  factor  at  each  mode, 
we  just  measure  bandwidth  of  the  half-power  points  of  the 
response  amplitude  from  the  detector  at  that  particular 
mode.  A  nondimensional  quantity,  quality  factor  Q,  can  also 
be  computed  from  knowledge  of  the  bandwith  by  the  relation! 

Q  =  -$r~  .  (4-4) 

4,  Experimental  Data 

In  order  to  verify  some  of  the  theoretical  predictions, 
we  have  kept  two  things  in  mind*,  the  experimental  condition 
under  investigation  and  the  physical  parameters  for  subhar¬ 
monics  to  be  excited*  The  experimental  data  are  recorded 
for  the  purpose  of  clarifying  how  good  the  theoretical 
analysis  is  for  a  real  physical  situation.  Those  data  are 
obtained  from  a  typical  interferometer  which  has  1.5  MHZ  and 
360  KHZ  X-cut  quartz  discs  for  its  boundary  walls, 

(a)  The  Resonance  Modes  of  the  Interferometer 
Table  4-1  lists  resonance  modes  of  an  interferometer 
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for  a  spacing  of  3*58  cm  at  26°C.  The  calculated  modes  are 
based  on  model  described  in  Chapter  III,  section  6,  and  com¬ 
puted  from  a  computer.  The  constants  for  the  characteristic 
equation  of  Eq.  (3-83)  used  during  computer  programming  ares 
R  =  10.2702 
fd=  1500.0 

f  =  360.0 
s 

f  =  21.0 
h 

The  measured  values  are  obtained  experimentally  from  the 
interferometer  used  in  subharmonic  observation.  The  proce¬ 
dure  of  such  data  taking  has  been  already  described  in  the 
previous  section. 

The  quality  factor,  Q,  associated  witn  the  resonance 
m^des  is  also  recorded  as  a  reference  to  indicate  the  sig¬ 
nification  of  such  a  particular  mode. 

(b)  Threshold  for  Subharmonic  Excitation 
There  are  two  parameters  for  setting  up  the  experiment 
for  observing  subharmoniess  one  is  the  length  of  the  inter¬ 
ferometer  and  the  other  is  the  driving  signal  frequency. 

The  dots  in  Fig,  4-iOa  are  the  data  for  the  subharmonic 
threshold  with  different  lengths  of  the  interferometer;  Fig, 
4-10b  is  a  record  showing  a  general  tendency  of  the  subhar¬ 
monic  threshold  to  change  as  a  function  of  the  driving  signal 
frequencies.  Those  data  are  taken  without  adjusting  other 
references  of  the  whole  experimental  system.  Since  the 


Pressure  (Bars)  Pressure  (Bars) 


1500  1520  1540  1560  1580  1600  1620 

Input  Signal  Frequency  (kHz) 

(b) 

FIG.  4-10  THRESHOLD  FOR  SUBHARMONiC  EXCITATION 
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subharmonic  frequency  components  change  with  interferometer 
length  and  driving  signal  frequency,  those  experimental 
points  do  not  have  the  same  frequency  components*  However, 
the  general  tendency  of  the  minimum  threshold  for  exciting 
the  subharmonics  can  be  indicated  by  the  solid  lines  in 
those  figures. 

(c)  Amplitude  and  Frequency  Components  of  Subharmonics 

Two  sets  of  experimental  conditions  have  been  carried 
out  in  a  detailed  study  of  the  amplitude  and  frequency  com¬ 
ponents  of  the  observed  subharmonics.  Fig.  4-lla  shows,  for 
a  case  of  subharmonic  of  one-half,  that  the  change  in  sub¬ 
harmonic  amplitude  measured  from  the  wave  analyzer  as  the 
driving  signal  intensity  increases.  Sometimes,  under  dif¬ 
ferent  conditions,  more  than  one  pair  of  subharmonics 
appear.  Fig.  4-llb  is  an  example  of  such  a  circumstance. 

The  frequencies  are  measured  by  an  electronic  counter  whose 
accuracy  is  about  *  5  HZ  in  our  present  measuring  range. 
Precision  of  the  data  on  the  reading  of  subharmonic  fre¬ 
quencies  is  dominated  by  the  stability  of  the  signal  gen¬ 
erator  which,  for  our  equipment,  is  about  +  10  HZ,  Hence, 
the  frequency  data  can  be  va3 id  up  to  +  15  HZ. 


armonic  Response  (mV) 


FIG.  4-11  INTENSITY  OF  SUBHARMONIC  RESPONSE 


5.  Discussion  of  the  Experimental  Results 

Some  physical  explanations  of  the  threshold  of  sub¬ 
harmonic  generation  are  gleaned  from  the  data  concerning 
changes  in  the  threshold  as  a  function  cf  length  of  the 
interferometer  as  shown  in  Fig0  4-10a.  When  the  distance 
between  the  boundary  walls  of  the  interferometer  is  small, 
the  loss  is  dominated  by  the  transducer  mounting  loss  which 
is  independent  of  the  interferometer^  length.  However, 
the  amount  of  medium,  which  provides  the  nonlinear  coupling 
mechanism,  is  proportional  to  the  space  between  the  boundary 
walls  of  the  interferometer.  These  two  factors  make  the 
threshold  for  exciting  the  subharmonics  decrease  with  the 
length  of  the  interferometer,  as  is  clearly  indicated  from 
the  experimental  data.  When  the  length  of  the  interfero¬ 
meter  becomes  very  large,  the  experimental  results  show 
that  the  subharmonic  generation  threshold  increases  with 
the  length  of  the  interferometer.  This  is  due  to  the  in¬ 
fluence  of  the  other  losses  such  as  the  viscosity  of  the 
medium,  beam  spreading,  etc.,  which  are  a  function  of  length 


The  fact  that  subharmonies  becomes  so  difficult  to  excite 
for  short  spacings  is  evidence  that  the  vibration  of  the 
boundary  walls  does  not  provide  an  important  contribution 
to  the  nonlinear  coupling  responsible  for  subharmonic  genera 
tion. 


In  Appendix  A,  we  have  carried  out  a  detailed  analysis 
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comparing  the  order  of  magnitude  of  the  nonlinear  coupling 
due  to  the  boundary  and  to  the  nonlinear  property  of  the 
medium. 


The  change  in  threshold  for  subharmonic  generation 
with  driving  signal  frequency,  as  indicated  in  Fig,  4-lOb, 
covers  a  quite  large  range  and  does  not  follow  a  sirrole 
relation.  This  is  due  to  the  reason  that  other  factors, 
such  as  dissipation  and  detuning  parameters,  should  also  be 
taken  into  consideration. 

Information  about  the  quality  factor,  G,  of  the 
interferometer  associated  with  resonance  modes  can  be  used 
to  calculate  the  threshold  for  subharmonic  generation  ac¬ 
cording  to  the  analysis  in  Chapter  Ilf.  From  Eq.  (3-?3). 
the  simplified  result  for  water  without  detuning  is: 


15  X  103 


(4-5) 


where  is  the  acoustic  pressure  in  bars,  G^,  and  Q2  are 
the  quality  factors  associated  with  subharmonic  modes  1  and  2 
respectively.  Based  on  the  data  of  Table  4-1,  the  calculated 
minimum  threshold  is  5*8  bars  which  is  close  to  what  we  have 
measured  in  Fig.  4-10. 

Though  the  interferometer  possesses  many  resonance 
modes,  in  general  only  a  few  pair  of  subharmonics  are  gener¬ 
ated  during  experiments.  This  can  be  explained  by  noting 
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that,  in  addition  to  the  threshold,  those  modes  also  have 
to  meet  a  matching  condition}  that  is,  the  sum  of  their  fre¬ 
quencies  should  also  equal  to  the  frequency  of  the  driving 
signal.  We  use  the  data  set  of  the  resonance  modes  from  Table 
4-1  to  check  such  a  situation.  Pig.  4-12  is  dotted  with 
points  corresponding  to  the  sum  of  frequencies  for  a  pair 
of  subharmonic  modes  in  the  range  from  1575  KHZ  to  1635  KHZ. 

Two  horizontal  lines  are  drawn  for  the  resonance  modes  at 
1595*3  and  1618,2  KHZ.  A  highly  intense  wave  can  be  estab¬ 
lished  if  we  drive  the  interferometer  at  these  frequencies. 

We  see  that  there  are  actually  only  few  points  near  these 
two  lines j  in  another  words,  not  many  subharmonic  pairs 
will  be  excited.  The  vertical  distance  between  those  points  to 
the  lire  gives  magnitude  of  the  resultant  detuning  and  numbers  in  pa¬ 
renthesis  are  the  geometric  means  of  the  quality  factors  associated 
with  those  particular  subharmonic  pairs.  Points  far  from 
the  line  may  be  excited  only  if  a  much  stronger  signal  is  pro¬ 
vided.  However,  the  final  amplitude  of  the  excited  subhar- 
raonic  will  depend,  not  only  on  the  quality  factor  for  that 
particular  mode,  but  also  on  its  detuning  from  the  actual 
resonance.  When  more  than  one  subharmonic  pair  is  excited, 
their  amplitudes  will  be  inversely  proportional  to  the  loss 
and  the  detuning.  The  experimental  results  shown  in  Pig. 

4-llb  give  such  evidence,  (Because  -f  ambient  temperature 
changes  due  to  the  different  times  for  data  collection,  there 


to  -  Sum  of  the  Mode  Number 
kN  =  Difference  of  the  Mode  Number 


4*12  SUBHARMONIC  PAIR  FREQUENCY 
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are  frequency  shifts  for  the  resonance  modes  in  Pig.  4-10 
from  the  original  modes  recorded  in  Table  4-1.  However,  we 
can  still  identify  them  by  the  general  direction  of  the 
shift.  The  corresponding  subharmonic  pairs  for  Figo  4-11 
are  labeled  with  s,  s^,  s2,  in  Pig.  4-12,  We  can  see  that, 
as  s^,  s2  both  have  high  Q*s  and  are  close  to  the  driving 
signal  frequency,  they  can  appear  at  the  same  time.  Since 
the  condition  for  measuring  the  resonance  frequencies  is  not 
the  same,  we  do  net  attempt  to  use  these  data  to  calculate 
the  threshold  and  subharmonic  intensity  for  the  fact  of 
detuning) . 

The  change  in  subharmonic  amplitude  with  driving  signal 

intensity  illustrates  a  surprising  discrepancy  from  results 

1 6 

reported  by  Bamberg  .  According  to  his  observations,  once 
the  subharmonic  is  excited,  its  amplitude  will  level  off 
regardless  of  ho w  intense  the  driving  signal  is  and  finally 
it  will  disappear  altogether  after  a  further  increase  in 
driving  signal  intensity »  We  feel  that  our  results  differ 
from  his  because  we  have  avoided  cavitation  during  our  ex¬ 
periments,  Since  gas  bubbles  consume  energy  through  the 
cavitation  process,  their  presence  will  certainly  prevent 
further  increase  in  the  amplitude  of  subharmonics.  This 
character  may  be  a  way  to  differentiate  between  subharmonic 
generation  with  cavitation  and  without  cavitation. 
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Chapter  V 
SUMMARY 


1.  Conclusions 

The  work  reported  here  has  been  an  attempt  to  develop 
a  mathematical  analysis  of  subharmonic  generation  in  acoustic 
systems  and  to  check  experimentally  the  adequacy  of  the  theo¬ 
retical  description  of  such  a  phenomenon. 

The  principal  conclusions  of  this  thesis  are  the 
following * 

(a)  Subharmonic  oscillation  is  a  manifestation  of  sys¬ 
tem  instability.  The  analysis  of  its  behavior  can  be  ob¬ 
tained  by  first  assuming  its  existence  and  then  determining 
the  conditions  under  which  i4-  occurs.  The  two-variable  (two- 
timing  )  perturbation  method  has  been  a  very  useful  mathe¬ 
matical  tool  in  carrying  out  this  study. 

(b)  In  a  conservative  system  of  three  oscillators 
with  a  single  nonlinear  element  whose  energy  function  is 
V=X1X2X^4,  if  the  sum  frequency  of  two  oscillators  is  close  to 
the  third  one,  there  exist  three  integral  constraints.  By 
using  these  constraints  in  the  relevant  four-dimensional 
phase  space  we  can  obtain  a  single  trajectory  to  describe 
the  system’s  motion, and  the  time  variable  can  be  found  by  an 
integral  along  that  trajectory. 

(c)  For  a  dissipative  system,  subharmonics  can  be 
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excited  by  a  sufficiently  strong  externa1  source,  since  the 
nonlinearity  of  the  system  nrovides  +.h<->  coupling  mechanism 
for  energy  conversion  to  the  lower  frequencies.  The  losses 
in  the  system  and  detuning  from  the  exact  frequency  matching 
condition  play  comparable  roles  in  determining  the  threshold 
required  for  subharmonic  excitation. 

(d)  Subharmonics  always  appear  in  pairs  except  in  the 
degenerate  case  of  the  subharmonic  mode  one-half.  lr.  a  dis¬ 
tributed  system,  such  as  the  acoustic  irterferometer,  mere 
than  one  pair  of  subharmonics  may  be  excited  provided  they 
meet  the  frequency  sum  rule  and  the  driving  source  is  in¬ 
tense  enough.  According  to  the  theory,  for  a  medium  with 
frequency  independent  losses,  the  subharmonic  of  one-half 
has  the  lowest  threshold  fer  excitation  for  a  given  detuning. 

(e)  Subharmonic  generation  has  been  observed  experi¬ 
mentally  in  an  acoustic  interferometer.  Because  cavitation 
was  carefully  avoided  during  the  investigation,  the  exoeri- 
mental  results,  such  as  the  threshold  level,  and  the  specific 
subharmonic  modes  excited,  are  consistent  with  the  theorexi- 
cal  predictions. 

2.  Future  Work 

The  theoretical  study  of  subharmonic  generation  can 
be  further  studied  for  nonlinear  element  with  other  form 
of  energy  functions.  For  example,  if  the  three  oscillator 


5-3 


model  possesses  a  sec  or.  d  energy  function,  then  there  will 
be  an  additional  phase  variable  and  one  of  the  three  origi¬ 
nal  conrtraints  dees  not  hold  anymore.  Therefore,  the  prob¬ 
lem  can  not  be  reduced  to  quadrature.  The  existence  of 
*  cher  types  of  constraints  remains  to  be  determined. 

The  one-dimensional  model  for  the  distributed  system 
should  be  extended  to  three-dimensions.  This  will  permit  a 
more  realistic  comparision  of  theory  and  experiment. 

Such  an  approach  will  also  allow  decomposition  cf  the  wave 
number  vectors  to  cover  the  more  general  case  when  the  prop¬ 
erties  of  the  medium  becomes  dispersive. 

On  the  experimental  front,  it  seems  that  tempera¬ 
ture  control  of  the  container  will  be  necessary  if  more 
precise  quantitative  measurements  are  desired.  A  0,2°  Centi¬ 
grade  change  will  cause  0.6  KHZ  shift  in  frequency  for 
modes  around  1.5  MHZ  -with  the  present  set-up. 

Liquids  other  than  water  can  also  be  used  for  observing 
subharmonic  generation.  The  experimental  data  will  supply 
the  information  on  the  nonlinearity  associated  with  differ¬ 
ent  liquids.  However,  purification  and  possible  cavitation 
deserves  investigation  before  considering  the  meaning  of  the 
experimental  measurements. 

It  may  be  interesting  to  extend  the  investigation  to 
subharmonic  generation  in  a  solid  medium  resonator.  With¬ 
out  the  complications  of  cavitation,  such  nonlinear  phe- 


i 

t  , 
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noseas  may  provide  a  different  practical  method  of  energy  conversion. 


Acknowledgments 

I  would  like  to  express  ay  appreciation  and  gratitude  to  Prof. 
R.  E.  Kronauer  for  his  guidance  throughout  this  research  project.  I 
wish  also  to  acknowledge  the  valuable  advice  from  Prof.  F.  V.  Hunt, 
the  helpful  suggestions  of  Dr.  R.  E.  Apfel,  the  patience  of  my  wife, 
Molly,  and  the  assistance  of  the  staff  of  the  Acoustics  Research 
Laboratory.  Finally,  I  would  like  to  acknowledge  the  financial 
support  of  this  research  by  the  Office  of  Naval  Research. 


A-I 


Appen*  Lx  A 

REFLECTION  FROM  A  V '  RATIONAL  dOUNDARY 


x*=|9 

p 

fa 


sintojt* 

k\  sinCwft'+k^x* ) 

^  —  ..0  It  X 


- »> 

R  *  sin(«o^t 3  -k  *  x  *+/) 


Fig.  A-l 


When  waves  reflected  from  a  vibrational  boundary,  their 
amplitudes  and  phases  will  change  according  to  the  boundary 
conditions.  Here  we  are  discussing  such  effects  and  make  a 
comparision  with  the  interaction  due  to  the  nonlinearity  of 
medium. 

Fig.  A-l  shows  that  a  wave  with  the  particle  velocity 
of  A£  sin  (co*t*  +  kjx*  )  is  incident  on  a  rigid  boundary 

whose  surface  moves  harmonically  with  a  displacement  as 
3^'  sin  (o0^t*)f  A  reflected  wave  can  be  found  by  matching 
the  boundary  conditions,  in  this  case,  it  is» 

R*  sin  (co£t*  -  k|x'  +  /’)  ,  (A-l) 

where  R'  =  -A£  , 

=  2k^’  sinoc^t*  . 

Asj  sin  Osina)  =  2J^O)sina  -  2J^O)sin3a  + 


ont 


nO)  -  l|)  -  jTfn-1  )i  C2*  4  *  *  *]  * 

t  A-  2 ) 

we  can  exnand  the  reflected  wave  by  considering  Th‘‘  t‘p+ 
order  of  only  (%’ is  a  small  quantity): 

R'  sin  (u?*t*  -  Kjx*  +  ff* ) 

-  -  A,  sin  (oj^t*  -  k£x*)  cos  (2^,kjsimo?  t9 ) 

-  AJ  cos  (oJit*  -  kjx*)  sin  , 2^*K ‘sinWJt ' ) 

=  -  A'  sin  («J*t*  -  K |x* )  [j0(2kj^')  *  2J2(2k^'  r.r .*2w?f  ' 

+  •*•••*•]  +  A{  cos  (o>*t»  -  k'x’)  [2J1t2k‘^*  )  *  c^‘ 

*  . 3 

«  -  Aj  sin  (w't'  -  k3*x*)  -  Aj  cos  (oo|t*  -  k{x') 

*  2k-j^‘sin  co^t ' 

=  -  AJ  -i»»  ^co^t *  -kjx')  -  kf^'A*  [sintw-'t1  -  u>jt*  A  <{v*  • 
+  sin  (w^t*  +  co£t*  -  kjx')J 

(A-3) 

The  amplitudes  of  new  frequency  components  (w^  -u>^)  &»'*i 
(to*  4  uxp  appear  ass 

Jfi- 

a  •-  lr  1  jb*  A  c  (AJ}  ) 

"2  "1*"1  ’  v"  ' 

and  in  the  nondimensional  notation  defined  in  Chapter  11!*  it 
becomes: 

Ag  =  .  (A-t) 

However,  the  vibrating  boundary  will  also  radi.ste 
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an  acoustic  v.*a/e  whose  particle  velocity  *isj 

U3  =  cos  c°3xl  *  (A-6) 

This  wave  v/ill  interact  with  the  reflected  wave  if  the 
c.cdui'i  possess;*®  a  nonlinear  property.  According  to  the 
in  Chapter  III,  a  new  frequency  component  of  oj^ 
or  -  oJ^  will  appear  v/ith  a  relation  ast 


(A-7) 

for  one  dimensional  problem, 

A^  =  6  ^(1+ P  , 

(A-3) 

v/here  L  is  the  distance  in  wavelengths  travelled  by  the  two 
interacting  waves. 

Comparing  Eq;s,  (A-?)  and  (A-5),  we  have: 


(A~9) 


Since-  and  are  in  same  order  and  P  is  about  6.5  for  v/ater, 
Eq,  (A-9)  indicates  that  is  L  times  larger  than  A £.  This 


is  proof  that  the  nonlinear  effect  due  to  medium  has 


more  influence  in  the  generating  new  frequency  components 
than  a  vibrational  boundary. 


Af.rrndix  B 


X 


k-h 


•*T/ATl*\.  Ht  2H2  Oi’ARAiiTERlSTJC  EQUATION’  FOR 
?rui*A  CL  uOCHS  I..  THE  THREE-SECT  I CV  CCWFOSITE  RESONATOR 


.  *.*  rcso-»n*i.  e  modes  cf  a  certain  resonator  ca-’  be  de- 
term  i~- a  frer  its  boundary  condition?*  However ,  for  the  one 
d'rer  anal  resonator  shown  in  Fir,  ?-2,  the  transmission 
5  inr  ;  .  T'.y  oa:.  ;  »>iy  for  solving  such  a  o^ob  1  en . 

formide^  rVr  section  of  tne  composite  resonator  as  a 
part  r{  ^  transmission  line.  Since  the  wall  at  rirht  side 
is  terminated  at  a  pressure  release  surface  (p  ~  0),  the 
acoustic  impedance  it  represents  to  the  left  will  be: 

x  -"v/'w  '  •  w  '  ^ 

v/here  all  the  notations  are  defined  as  before. 

The  middle  section  (medium)  will  transfer  Z^  to  its 
left  side  with  a  new  impedance  value  as: 


Z  i  c '  tan  k '  S 

z  =  &  *  c  *  — _ _ _ _ _ 

2  ^  c  ^'o'  +  tan  K'S 


(B-2) 


The  resultant  impedance  appearing  at  the  left,  wall  will 
then  become: 


z  -r.r.  Z2  4  ipw  tan  kJg 

^  '  tat  *  X 


(B-3) 


w  w  p’c’  ,iZ0  tan  I'D  * 
v  w  w  2  w 

As  tho  left  wall  also  terminates  v/ith  a  pressure  ro- 


lor.s*'  surface. 


at  resonance 


we  should  have: 


A  . 


t 


Z3  =  0  . 

With  the  relations  of  Eqs.  (B-t)r  ( 
we  gets 

[R(tan  k^D  +  tan  k^H)  +  tan  k*S  -  R2  tan 
/(R  -  R2  tan  k'H  tan  k*S  -  R  tan  k*K  tan 

T9  W 

xian  k£P) 

»  0  , 

or  simply, 

R(tan  k*D  +  tar  k*K)  +  tan  k*S  -  R2  tan 
Vt  w 

xtan  k*S 

=  0  , 


A- 

(B-4) 

H-?}.  (F->3). 

k,*D  tan  k*K  tan  k* 
w  w 

k*D  -  tan  k*S 

wi 

(B-5) 

k^D  tan  k£H 

<B-5> 


which  is  just  the  characteristic  equation  0-63)* 


Anoonriix  C 
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Table  C-l  lists  the  computer  program  for  carrying  out 
the  actual  computation  process  of  Eq.  (C-3) ♦  As  f_  5s  much 
smaller  than  and  f^t  the  separation  between  two  succes¬ 
sive  roots  will  not  exceed  f_.  We  use  0*65^  for  the  ini- 
tial  value  during  iteration. 

In  this  program,  the  region  near  singular  points  of 
the  tangent  functions  is  avoided  to  prevent  the  computer  from 
overflowing.  Any  missed  roots  are  determined  by  interpolating 
from  their  neighboring  roots. 

Notations  used  in  the  program  follow  closely  to  what 
are  defined  in  Pig.  3-2.  The  final  output  solution  Y(I)  is 
the  frequency  of  the  resonance  mode  in  units  of  KKZ,  and  G(I) 
is  the  actual  value  of  the  function  F{f)  when  f  =Y ( I ) , 
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